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Chapter  1 
Introduction 


The  main  objectives  of  this  study  were  to  develop  a  re-usable  flight  control  law,  advance 
control  allocation  technology  and  apply  the  control  law  to  a  challenging  configuration.  A 
modular  flight  control  architecture  is  designed  that  will  maximize  software  reusability  for 
future  control  law  development  efforts  and  allow  efficient  advanced  technology  insertion  as 
it  becomes  available.  A  control  allocation  method  is  developed  that  optimally  utilizes  all 
available  control  power  without  violating  actuator  limits.  A  tailless  fighter  configuration  was 
chosen  for  application  due  to  the  operational  benefits  and  flight  control  law  development 
challenges  inherent  in  a  tailless  configuration  [6] . 

A  modular  flight  control  system  design  approach  is  used  so  that  promising  technologies 
from  feedback  control  synthesis  and  control  allocation  may  be  evaluated  interchangeably 
without  complete  control  system  re-design.  The  control  law  is  based  upon  feedback  lin¬ 
earization  or  dynamic  inversion  control  theory  and  follows  closely  with  the  MACH  control 
philosophy  in  the  MCT  guidelines  [25]. 

Control  allocation  algorithms  mtist  efficiently  use  nonlinear  multi-axes  control  effectors 
to  achieve  current  maneuverability  requirements  of  fighter  aircraft.  The  control  allocation 
method  should  exploit  control  effector  redundancy  to  optimize  mission  segment  objectives 
such  as  minimize  drag  or  wing  loads.  The  control  allocation  must  also  optimize  large  com¬ 
mand  response  and  prevent  departures.  Although  not  specifically  addressed  in  this  study, 
low  computational  requirements  is  a  very  important  characteristic  of  the  control  allocation 
method.  A  tradeoff  exists  between  control  allocation  computational  efiiciency  and  perfor¬ 
mance,  and  tools  exist  to  analyze  this  tradeoff  [24].  Computational  efficiency  provisions  have 
been  added  to  existing  methods  [16]  at  the  expense  of  performance.  The  control  allocation 
method  developed  in  this  study  is  a  refined  version  a  linear  program  based  method  [10]  which 
optimally  exploits  effector  redundancy,  optimizes  large  command  aircraft  response  and  uses 
all  the  available  control  power. 

Tailless  aircraft  configurations  have  gained  recent  support  due  to  inherent  increase  in 
stealth  and  decreases  in  weight  and  drag.  Flight  control  challenges  for  this  configuration 
include  multi-axes  instabilities  and  nonlinear  control  effectors.  The  operational  advantages 
and  the  flight  control  challenges  drove  the  current  study  toward  a  tailless  fighter  configuration 
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This  docTiment  is  organized  as  follows.  Chapter  2  describes  background  material  regard¬ 
ing  the  aircraft  configuration  and  dynamic  inversion  control  theory  as  it  applies  to  flight 
control.  Chapter  3  contains  a  modular  flight  control  law  that  is  applied  to  the  tailless  air¬ 
craft  configuration  in  Chapter  2.  Chapter  4  provides  a  detailed  analysis  of  the  tailless  aircraft 
flight  control  law.  Chapter  5  summarizes  the  most  important  results  and  conclusions  of  the 
flight  control  development  and  analysis. 
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Chapter  2 
Background 


2.1  ICE  Configuration  101-3 

The  vehicle  in  this  study  is  a  tailless  fighter  configuration  developed  under  the  Innovative 
Control  Effectors  (ICE)  program  [15].  Although  there  were  multiple  configurations  studied 
in  the  ICE  program,  only  configuration  101-3  was  chosen  in  the  current  study  to  evaluate 
an  advanced  control  law.  Configuration  101-3  was  the  most  promising  from  an  aerodynamic 
perspective  while  being  the  most  challenging  from  a  control  law  development  perspective. 
The  ICE  configuration  101-3  is  shown  in  Figure  2.1. 

2.1.1  Control  Effectors 

The  control  effectors  for  configuration  10T3  include  elevens,  pitch  flaps,  all  moving  tips, 
thrust  vectoring,  spoiler  slot  deflectors,  and  outboard  leading  edge  flaps.  The  conventional 
control  effectors  are  defined  as  the  elevens,  pitch  flap  and  leading  edge  flaps.  The  innovative 
control  effectors  are  defined  as  the  thrust  vectoring,  all  moving  tips  and  spoiler  slot  deflectors. 
Pitch  and  yaw  thrust  vector  coupling  results  from  the  non-rectangular  nature  of  the  thrust 
vector  limits.  The  challenges  of  the  all  moving  tips  and  spoiler  slot  deflectors  are  zero  lower 
deflection  limits,  strong  multi-axes  effects  and  effector  interactions. 

2.1.2  Simulation  Model 

A  generic  simulation  environment  is  used  for  ICE  configuration  simulation  model  develop¬ 
ment.  The  simiilation  model  contains  ordinary  differential  equations  of  the  aircraft  dynamics 
of  the  following  form 

X  =  f{X,u^)  (2.1) 

where  X  are  the  aircraft  states,  including  actuator  states,  and  n*'  are  the  control  effector 
commands.  The  right  hand  side  of  the  differential  equations  {f{X,  u))  is  computed  using  the 
aerodynamic  model  developed  from  Phase  I  wind  tunnel  experiments,  a  static  engine  model, 
actuator  models  and  a  standard  atmosphere  model.  Although  sensor  noise,  turbulence  and 
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Figure  2.1;  ICE  Configuration  101-3 
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gusts  are  also  modeled,  no  results  are  presented  since  only  a  cursory  analysis  was  completed 
with  these  effects. 

The  ICE  program  [15]  developed  hinge  moment  data  which  was  used  to  size  the  actuators. 
However,  simulation  analysis  showed  that  maximum  hinge  moments  were  exceeded  during 
high  speed  maneuvering.  The  all  moving  tip  pivot  axis  was  initially  located  such  that  effector 
hard-overs  resulted  if  maximum  hinge  moments  were  exceeded.  The  hinge  moment  effects 
were  neglected  in  this  study  due  to  the  need  of  a  more  rigorous  actuator  sizing  and  all 
moving  tip  pivot  location  analysis.  The  simulation  model  of  each  control  effector  actuator 
included  fourth  order  dynamical  models  with  deflection  and  rate  limits.  The  properties  of 
each  actuator  are  given  in  Table  2.1. 

The  thrust  vector  nozzle  deflection  limits  are  conical  about  the  zero  deflection  engine 
centerline.  Note  the  low  phase  crossover  frequency  of  the  thrust  vector  nozzle  actuator. 

2.2  Dynamic  Inversion  Flight  Control  Theory 

The  aircraft  rigid  body  equations  of  motion  in  eq.(2.1)  may  be  approximated  by  ordinary 
differential  eqTiations  that  are  affine  in  the  control  if  actuator  dynamics  and  control  effector 
interactions  are  neglected 

X  =  a{x)  +  b{x)u  (2.2) 

where  .r  €  IR”"  are  the  states  (excluding  actuator  states),  u  G  IR"*  are  the  control  effector 
deflections  (or  actiiator  states)  and  a(.r)  and  columns  of  b{x)  are  smooth  vector  fields  on  IR”. 
All  rigid  body  states  are  typically  measured  and  available  for  feedback  control. 

A  set  of  feedback  control  command  variables  is  defined  as  a  generally  nonlinear  function 
of  the  states 

y  =  c(x)  (2.3) 

where  y  G  IR”"  is  the  vector  of  command  variables.  It  is  ass\imed  that  there  are  at  least  as 

many  controls  as  command  variables,  i.e.  m  >  Uy.  There  are  typically  at  least  3  controls  for 

the  3  command  variables  for  manual  flight  control.  It  is  assumed  that  the  command  variables 
in  eq.(2.3)  for  the  system  in  eq.(2.2)  have  a  well-defined  vector  relative  degree  with  unity 
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elements  {ri,  •  -  • , r„ J  =  {1,  •  •  • ,  1}.  Typically  each  manual  flight  control  command  variable 
has  unity  relative  degree,  however  the  theory  is  easily  extended  for  other  than  unity  relative 
degree.  Given  these  assumptions,  there  always  exists  a  diffeomorphic  state  transformation 
to  the  following  convenient  form 


z 

y 


az{z,y) 

0’y{^,y) 


+ 


{z,y) 

<y{z,y) 


(2.4) 


where  2  are  the  uncommanded  states.  A  well-defined  vector  relative  degree  implies  that  time 
derivatives  of  all  command  variables  are  affected  by  at  least  one  control  and  that  by  has  rank 
equal  to  Uy.  It  is  further  assumed  that  each  control  affects  at  least  one  command  variable 
time  derivative  which  is  typical  for  manual  flight  control. 

The  feedback  linearization  or  dynamic  inversion  control  law  has  two  main  parts  listed 
below 


U  =  p{dy) 

dy  =  -ay{z,y)-^v{y,f^).  (2.5) 

The  first  part  is  the  control  allocation  (/?)  which  is  a  mapping  from  generalized  controls  {dy) 
to  actual  controls  (u)  such  that  the  following  relationship  holds 

by{z,y)p{dy)  =  dy  .  (2.6) 

The  second  part  is  the  generalized  control  law.  This  part  consists  of  deaugmentation  of 
the  natural  command  variable  dynamics  (-ay(z,y))  and  augmentation  of  some  desired  com¬ 
mand  variable  dynamics  (v(y,y^‘‘))  that  depend  on  the  control  system  requirements.  The 
reference  command  vector  is  y^'^. 

The  nominal  closed  loop  is  defined  by  the  following  system  with  a  cascade  structure 


2 

fz{z,y) 

.  y . 

v{y,y<^^)  J 

fz(z,y)  =  a^(z,y)  +  b,(z,y)p(-o.y(z,y)-{-v(y,y^‘‘))  (2.7) 

It  is  may  be  shown  [26]  that  this  system  is  asymptotically  stable  if  the  command  variable  sub¬ 
system  (y)  is  asymptotically  stable  and  the  complementary  subsystem  (z)  is  asymptotically 
stable  on  the  zero  command  variable  manifold  which  is  equivalent  to  asymptotic  stability  of 
the  following  system 

i  =  0.2(2, 0)  -f-  6^(2, 0)tt(2)  (2.8) 

where  u(z)  is  the  control  that  maintains  zero  command  variables 

u(z)  =  {u\ay{z,  0)  -t-  by{z,  0)u  =  0}  .  (2.9) 

Note  that  the  system  in  eq.(2.8)  is  often  referred  to  as  the  zero  dynamics,  and  asymptotic 
stability  of  the  zero  dynamics  is  equivalent  to  the  minimum  phase  property.  Further  note 
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that  asymptotic  stability  of  the  command  variable  dynamics  is  guaranteed  by  choosing  a 

stabilizing  j  •  ui  •  (o 

Stability  of  the  zero  dynamics  depends  on  the  choice  of  command  variables  m  eqd^-^j- 

There  are  methods  to  choose  the  command  variables  to  optimize  stability  robustness  of  the 
zero  dynamics  [27].  However,  the  choice  of  the  command  variables  also  affect  disturbance 
rejection  and  robustness.  A  command  variable  choice  that  optimally  rejects  disturbances  in 
the  presence  of  uncertainties  may  result  in  unstable  zero  dynamics.  In  flight  control,  the 
zero  dynamics  are  typically  stable  for  an  optimal  robust  disturbance  rejection  choice  of  the 

command  variables. 
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Chapter  3 

Flight  Control  Structure 


The  proposed  flight  controller  (MODCON)  for  the  ICE  aircraft  in  Section  2.1  is  shown  in 
Fig.  3.1. 


measurements 


Figure  3.1:  Controller  Structure 

It  has  a  modular  structure  that  encourages  insertion  of  independent  technology  advances 
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in  feedback  control,  modeling  and  control  allocation.  The  aircraft  MODEL  module  provides 
vehicle  specific  data,  such  as  stability  and  control  derivatives,  and  requirements  data,  such 
as  desired  bandwidths,  to  all  other  control  law  modules.  The  feedback  and  feedforward 
aiigmentation  (FFA)  module  generates  control  law  commands  that  provide  desired  stability 
and  flying  qualities  properties  in  a  manner  that  is  robust  to  uncertainties.  The  control 
allocation  (CONALL)  module  computes  actuator  commands  that  optimize  a  mission  segment 
objective  for  achievable  control  law  commands  and  optimally  limits  unachievable  control 
law  commands.  The  command  shaping  (CMDSHP)  module  provides  proper  feedforward 
characteristics  such  as  stick  force  gradients,  command  limiters  and  command  filters. 


3.1  Command  Shaping 

The  command  shaping  module  provides  appropriate  command  augmentation  properties.  The 
F-16  was  used  as  a  requirements  baseline.  Therefore,  a  similar  side  stick  shaping  structure 
to  the  F-16  was  used  for  the  ICE  aircraft.  The  command  shaping  module  consists  of  three 
sub-modules:  stick  gradient,  command  limiter  and  command  filter.  The  stick  force  gradient 
maps  force  inputs  in  pounds  to  CV  commands  in  degrees  per  second.  The  command  limiter 
restricts  the  CV  commands  to  help  prevent  departure.  The  command  Alter  provides  desired 
transient  command  response. 


3.1.1  Stick  Gradient 

The  stick  gradient  provides  desirable  mappings  between  pilot  stick/pedal  forces  and  flight 
control  commands.  The  piece-wise  linear  structure  for  the  stick  and  pedal  gradients  is  shown 
in  Fig.(3.2).  The  gradient  parameters  for  each  axis  are  supplied  by  the  on-board  model  and 
are  given  in  Section  3.4. 

3.1.2  Command  Limiter 

Only  a  cursory  design  of  the  command  limiter  was  completed  due  to  time  limitations.  The 
command  limits  were  chosen  to  be  of  the  following  form 

CVmax  =  min  {kiq,  k^)  (3-1) 

where  q  is  the  dynamic  pressure.  The  command  limit  parameters  {ki,k2)  for  each  axis  are 
supplied  by  the  on-board  model  and  are  given  in  Section  3.4. 


3.1.3  Command  Filter 

The  command  filter  in  each  control  channel  has  the  following  first-order  form 


(3.2) 
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Command 


Figure  3.2:  Stick  Gradient  Structure 

where  Uc  is  the  desired  feedback  control  bandwidth  and  Ufq  is  the  desired  flying  qualities 
command  bandwidth.  The  command  bandwidth  Ufq  allows  the  capability  of  fast  command 
response  without  sacrificing  closed  loop  stability  margins.  The  command  filter  parameters 
for  each  axis  are  supplied  by  the  on-board  model  and  are  given  in  Section  3.4. 


3.2  Feedback  and  Feedforward  Augmentation 

Dynamic  inversion  synthesis  [10,  22]  is  used  for  feedback  and  feedforward  augmentation. 
Dynamic  inversion  consists  of  deaugmentation  of  the  natural  dynamics  of  a  specified  set 
of  command  variables  followed  by  robust  feedback  stability  augmentation  and  feedforward 
command  variable  augmentation.  The  generalized  dynamic  inversion  control  law  is  developed 
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in  Section  2.2  and  is  given  by 


^des  ^ 

v{y,y^'')  =  ujcXi-uJcy  +  ujJcy^^^^ 

X,  =  -U,hy  +  ujcfiV^  (3.3) 

where  y  are  the  command  variables,  are  the  generalized  control  commands,  —o.y(x)  is 
the  deaugmentation  part  and  ?>(?/,  y"”'*)  represents  the  augmentation  part  which  consists  of 
proportional-integral  stability  augmentation  and  command  augmentation.  Dynamic  inver¬ 
sion  inverts  the  command  variable  dynamics  and  proportional/integral  feedback  generates 
a  desired  robust  loop  shape.  Proportional  gains  are  used  in  the  feedforward  path  to  obtain 
desirable  flying  qualities.  Note  that  is  the  reference  command  that  is  to  be  tracked  by 
the  command  variables.  The  augmentation  parameters  are  supplied  bj^  the  on-board  model 
and  are  given  in  Section  3.4. 


3.3  Control  Allocation 

By  design,  the  control  law  in  eq.(3.3)  provides  desired  loop  and  command  response  properties 
in  terms  of  however  the  reduced  dimension  controls  {dy^^)  must  still  be  resolved  into 
actual  control  effector  commands  The  control  allocation  function  in  eq.(2.5)  provides 
the  mapping  from  reduced  dimension  controls  to  actual  control  commands,  u'^  =  p[dy^^).  The 
control  allocation  module  computes  actuator  commands  that  achieve  the  desired  control  law 
command  if  possible,  or  a  limited  control  law  command  if  necessary,  to  satisfy  actuator 
limits. 

The  control  allocation  problem  has  been  stated  many  times  in  various  forms  [18,  30,  11, 
25,  10,  20].  The  unlimited  control  allocation  problem  is  to  find  u  such  that  ByU  -  for 
all  df\  For  redundant  control  effector  suites,  there  are  many  sohitions  to  this  unlimited 
control  allocation  problem.  A  unique  solution  is  found  by  introducing  a  performance  index 
that  forces  the  redundant  effectors  to  satisfy  additional  objectives  such  as  minimal  radar 
cross  section,  drag,  or  wing  loads.  However,  since  there  are  limits  on  the  control  effectors, 
not  all  df^  are  achievable.  Therefore,  df^  may  need  to  be  clipped  or  limited  such  that 
u  =  p{dy^^)  does  not  violate  actuator  limits. 

Past  control  allocation  research  includes  algorithms  that  minimize  control  deflections 
[11]  or  drag  [17]  and  limit  unachievable  control  law  commands  by  preserving  its  direction 
[1,  18,  4].  However,  it  has  been  shown  that  preserving  the  direction  of  the  command  may 
unnecessarily  degrade  flying  qualities  [10]. 

This  section  poses  the  control  allocation  problem  as  a  nonlinear  constrained  optimization 
problem.  Simplifying  assumptions  are  stated,  and  the  control  allocation  problem  is  trans¬ 
formed  into  a  standard  linear  program.  An  algorithm  is  designed  to  solve  the  constrained 
control  allocation  problem  using  a  linear  program  solver  as  its  basis. 


Air  Force  Research  Laboratory 


12 


3.3.1  Single-branch  Control  Allocation 

The  control  allocation  problem  with  no  actuator  limits,  or  unlimited  control  allocation  prob¬ 
lem,  may  be  cast  into  the  following  optimization  problem  with  a  nonlinear  objective  and 
linear  equality  constraints  that  enforce  the  control  law  command 

nun  J  =  /(u)  subject  to  ByU  =  .  (3.4) 


Analytical  solutions  exist  to  the  problem  in  eq.(3.4)  for  certain  forms  of  the  objective.  For 
example,  a  quadratic  objective  results  in  the  common  generalized  right  inverse 


J  =  ^  -  Upref)^  W  {U-  Upref) 

=  W-'B^  (B,W-'Bjy'  +u^.,  (3.5) 

Limits  on  the  actuators  complicate  the  solution  of  the  control  allocation  problem.  Nonlinear 
inequality  constraints  are  added  to  prohibit  violation  of  actuator  limits 


mm 

u 


J  =  f{u) 
subject  to 


ByU  = 

u  <  u  <  u 


(3.6) 


Rate  and  position  limits  are  accounted  for  by  defining  u  and  u  as  the  most  restrictive  actuator 
constraints 


u  =  min(u„,,  ATur) 

u  =  max(u(,  —  ATur)  (^•'^) 

where  is  the  actuator  tipper  position  limit  vector,  ui  is  the  actuator  lower  position  limit 
vector,  Ur  is  the  no-load  actuator  rate  limit  vector  and  AT  is  the  digital  flight  control  system 
update  rate. 

As  stated,  the  problem  in  eq.(3.6)  may  be  infeasible  since  there  is  no  guarantee  that  the 
equality  constraints  can  be  satisfied  in  the  presence  of  the  inequality  constraints.  Additional 
constraints  may  be  added  to  insure  feasibility 

mm  j  =  j{u)-\  (3.8) 

r  ByU  = 

subject  to  <  d  <  A  <  1 
[  u<u<u 

assuming  that  “u  =  0  satisfies  the  inequality  constraints.  This  assumption  is  not  restrictive 
if  increments  from  past  values  are  used  for  all  variables.  Increments  are  natural  for  digital 
control  system  implementation. 

The  scalar  A  allows  the  algorithm  to  relax  the  equality  constraint  if  df^  is  not  achievable 
and  therefore  provides  an  indication  whether  is  achievable.  If  A  =  1  then  is  achiev¬ 
able,  otherwise  it  is  not.  In  effect  A  limits  d^®®,  if  necessary,  by  reducing  its  magnitude  and 
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maintaining  its  direction.  This  optimization  formnlation  is  control  law  command  feasibility, 
control  deficiency  and  control  sufficiency  in  a  single  step.  Control  law  command  feasibility  is 
indicated  by  A  =  1  and  infeasibilit.y  is  indicated  by  A  1 .  Control  sufficiency  optimization  is 
achieved  through  an  appropriate  choice  of  f{u).  Control  deficiency  is  achieved  by  maximiz¬ 
ing  the  magnitude  of  the  control  law  command  A  while  maintaining  direction.  However,  it 
is  has  been  shown  that  preserving  the  direction  of  the  command  may  unnecessarily  degrade 
performance  and  thus  is  not  optimal  [10].  A  multi-branch  formulation  allows  more  fiexibility 
and  optimality  during  control  deficiency. 

3.3.2  Multi-branch  Control  Allocation 

Redundant  control  effectors  and  actuator  limits  require  essentially  two  parts  for  every  control 
allocation  algorithm:  control  sufficiency  to  provide  solution  uniqueness  in  the  presence  of 
control  redundancy  and  control  deficiency  to  prohibit  violation  of  actuator  limits.  Sufficiency 
or  deficiency  is  determined  through  a  third  part,  control  law  command  feasibility.  Although 
the  previous  section  shows  that  single-branch  algorithms  may  be  constructed,  a  multi-branch 
formulation  allows  more  fiexibility  and  optimality  during  control  deficiency. 


Control  Law  Command  Feasibility  and  Control  Deficiency 

The  control  law  command  feasibility  and  deficiency  functions  are  still  achieved  in  a  single 
step.  The  first  step  in  the  multi-branch  control  allocation  algorithm  is  determining  if  the 
control  law  command  is  feasible.  If  it  is  feasible,  control  sufficiency  provisions  are  made 
in  a  separate  control  allocation  branch  and  if  it  is  deficient,  control  deficient  provisions  are 
automatic.  Control  law  command  feasibility  may  be  determined  by  solving  the  following 
nonlinear  optimization  problem 

mm  J  =\\ByU-df^\\^  (3.9) 

subject  to  u<u  <u 


Note  that  the  1-norm  is  chosen  to  be  amenable  with  linear  program  solvers.  The  feasibility 
optimization  may  be  transformed  into  a  standard  linear  program 


mm  J  = [ 0  •  •  •  0  1 

Us 

—u 

subject  to  w 

-ByU  -f  Us 
ByU  Us 


(3.10) 


where  Us  is  a  slack  variable.  If  J  =  0,  the  control  law  command  is  feasible  otherwise  it 
is  infeasible.  Linear  programming  is  chosen  due  to  the  ubiquitous  solvers  commercially 
available. 
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The  solution  to  the  problem  posed  in  eq.(3.10)  may  also  be  used  during  control  deficiency. 
Since  the  objective  of  this  optimization  is  to  minimize  the  slack  variables,  Ug,  the  control 
deflection  solution,  u*,  provides  an  achieved  dy  that  is  closest  to  d*"'  in  the  1-norm  sense. 
Weights  may  be  added  for  greater  flexibility  in  achieving  optimality  during  control  deficiency 

mm  j^wjug  (3.11) 

subject  to 

Wj  e  IR”^' 

Ug  e  IR"*' 

Summarizing,  solution  to  eq.(3.11)  with  any  linear  program  solver  may  be  used  to  de¬ 
termine  feasibility.  If  J  =  0,  the  control  law  command  is  feasible  and  control  sufficiency 
provisons  must  be  made  which  is  the  subject  of  the  next  section.  If  J  7^  0,  the  control  com¬ 
mand  is  infeasible.  Further,  the  solution  from  eq.(3.11)  provides  optimal  control  deficiency 
in  the  sense  of  minimizing  the  weighted  1-norm  distance  between  ByU*  and  dy^^  where  u*  is 
the  optimization  solution. 

Control  Sufficiency 

If  the  solution  of  eq.(3.11)results  in  J  =  0,  there  is  sufficient,  and  possibly  excess,  control 
power  to  achieve  the  control  law  command,  df^.  The  excess  control  power  is  used  to  op¬ 
timize  mission  segment  objectives  such  as  minimum  drag.  The  following  linear  program  is 
formulated  for  control  sufficiency 

min  J  =  WyUg 

Ug 

—u 

subject  to  u 

-u  +  Ug 

U  +  Ug 

ByU  =  df  ^ 

Wj  €  IR”* 

Us  elR'” 

where  Upre/  and  Wu  are  chosen  depending  on  the  mission  segment  objective. 

Six  control  allocation  modes  were  implemented  to  address  specific  mission  objectives.  The 
minimum  control  deflection  control  allocation  mode  is  designed  to  minimize  deflections  of  all 
the  control  effectors  which  is  a  good  approximation  of  multiple  objectives  such  as  minimum 
drag  and  minimum  actuator  power.  The  minimum  drag  mode  is  designed  to  reduce  drag  by 
adding  a  higher  fidelity  on-board  model  of  drag  due  to  control  deflections.  The  minimum 


(3.12) 


■  0 
—u 
>  u 

'^pref 

Upref 
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wing  loading  mode  is  designed  to  reduce  wing  bending  loads  by  minimizing  the  use  of  the  all 
moving  tips,  the  most  outboard  effectors.  The  minimum  radar  signature  mode  is  designed 
to  reduce  radar  signature  by  minimizing  the  iise  of  spoiler  slot  deflectors  which  qualitatively 
give  higher  radar  signatures  than  the  other  effectors.  The  minimum  thrust  vector  mode  was 
designed  to  emphasize  the  use  of  the  aerodynamic  effectors.  The  final  mode  is  the  null  space 
injection  mode  which  is  described  in  the  next  section. 

Null  Space  Injection  Mode  The  eventual  goal  is  to  integrate  the  linear  program  control 
allocation  with  on-line  adaptation  to  provide  robustness  to  large  errors  and  optimal  per¬ 
formance  [8].  A  fundamental  requirement  for  adaptive  control  is  persistence  of  excitation 
and  a  well-conditioned  regressor  matrix  [13].  Therefore,  activity  of  all  control  effectors  is 
necessary  to  accurately  estimate  the  control  effectiveness  on-line.  Furthermore,  it  has  been 
experimentally  observed,  that  columns  of  the  regressor  matrix  should  be  of  the  same  order 
of  magnitude  at  each  point  in  time,  not  just  in  the  I2  norm  over  the  window.  This  is  a  noise 
effect,  viz.,  small  deflections  yield  low  SNR,  which  is  responsible  for  inaccurate  estimates. 
Simulation  analysis  showed  that  the  linear  program  control  allocation  typically  results  in 
most  effectors  following  their  preferences  while  using  only  a  subset  of  the  control  effectors  to 
achieve  the  control  law  command.  This  characteristic  would  make  identification  of  control 
effectivenesses  of  each  surface  impossible  since  the  deflections  may  be  correlated  and/or  may 
not  be  of  the  same  order  of  magnitude  for  all  preferences  (wpre/)  and  associated  weights 
(VF,,).  For  example,  static  preferences  and  weights  (•Upre/=constant,  VF„=constant)  result 
in  highly  correlated  solutions  and  high  variability  in  magnitude.  Therefore,  in  this  section 
alternative  control  preferences  and  weights  are  tised  to  assist  in  estimating  accurate  control 
derivatives.  The  motivation  for  the  injection  control  allocation  mode  is  to  maximize  the 
condition  number  of  a  regressor  matrix  consisting  of  control  deflections. 

Actuator  commands  of  uniform  distribiition  are  achieved  with  the  following  weighted 
pseudo- inverse  control  preference 

Ujrref  =  W-^  B];  {ByW-^  B^)-^  .  (3.13) 

Evidently,  these  commands  will  be  highly  correlated.  Decorrelation  is  accomplished  through 
random  variations  of  the  weighting  matrix  in  eq.(3.13) 

W  =  WWr 

Wr  =  diag(10""bl0’'"^’-)  (3.14) 

where  IF  is  a  nominal  diagonal  weighting  matrix  that  scales  each  effector  based  upon  de¬ 
flection  limits  to  equally  distribute  commands,  and  rv  is  a  vector  of  uniformly  distributed 
random  variables  between  -1  and  1. 

Most  control  derivative  identification  enhancement  may  be  accomplished  using  the  con¬ 
trol  preference  in  eq.(3.13)  and  random  weighting  in  eq.(3.14).  Although  Uprej  in  eq.(3.13) 
with  weight  W  in  eq.(3.14)  is  completely  decorrelated  and  equally  distributed,  it  is  just  a 
preference.  Unequally  distributed  actuator  commands  may  still  be  prevalent  due  to  the  na¬ 
ture  of  extremals  in  linear  programming.  So  stronger  linear  program  solution  enforcement 
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to  Ujrref  will  HO  doTibt  deter  unequal  distribution  and  correlation.  This  is  accomplished  using 
the  following  dynamic  control  weight  in  the  linear  program  in  eq.(3.12) 


Wu  =  diag( 


|uOi  -  nOpi  I  |n02  -  uOj^  | 


(3.15) 


where  uO  is  the  effector  past  value  vector  and  uOp  is  the  effector  preference  past  value  vector. 
Essentially,  this  dynamic  weighting  drives  the  effectors  farthest  from  their  preference  hardest 
toward  their  preference. 


Summary 

The  multi-branch,  linear  program-based  control  allocation  method  is  summarized  as  follows. 
A  By,  ui,  Uu,  Ur,  Wd,  Wy  and  Upref  is  supplied  by  the  MODEL  module  and  a  is  supplied 
by  the  FFA  module.  The  linear  program  in  eq.(3.11)  is  solved  using  these  data  {By,  ui,  w„, 
Ur,  dy^^,  Wd)  giving  the  optimal  solution  u*.  If  J{u*)  ^  0,  then  u^  —  u*  gives  the  closest 
moment  ByU'^  to  in  the  sense  of  the  performance  index  in  eq.(3.11)  without  violating  the 
actuator  limits.  If  J  =  0,  the  linear  program  in  eq.(3.12)  is  solved  using  these  data  {By,  ui, 
Uy,  Ur,  df\  Upref-,  ^u)  giving  u* .  Setting  u'"  =  u*  gives  the  desired  moment  ByU*  =  df^  that 
achieves  the  mission  segment  objective  in  the  sense  of  the  performance  index  in  eq.(3.12) 
without  violating  the  actuator  limits. 


3.4  On-board  Model 

The  purpose  of  the  on-board  model  control  law  module  is  to  compute  all  aircraft  configu¬ 
ration  and  control  law  requirements  data.  For  inner  loop  flight  control,  the  most  notable 
aircraft  configuration  data  are  the  stability  and  control  derivatives.  The  most  notable  con¬ 
trol  law  req\urements  data  are  desired  flying  qualities  parameters  and  feedback  bandwidths. 
The  aircraft  configuration  model  caused  the  on-board  model  to  be  the  control  law  module 
that  required  the  most  development  effort,  and  it  is  the  module  that  still  needs  the  most 
refinement,  especially  at  high  speeds. 

3.4.1  Requirements  Model 

The  requirements  model  contains  all  functions  and  parameters  that  pertain  to  control  law 
requirements.  Data  include  nominal  stability,  nominal  performance,  robustness  and  opti¬ 
mization  requirements. 

Command  Variables  Requirements 

A  set  of  feedback  command  variables  are  defined  such  that  acceptable  levels  of  robustness  to 
uncertainty  and  disturbance  rejection  are  achieved.  The  following  command  variables  were 
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chosen  based  on  control  law  efforts  for  a  similar  class  of  fighter  aircraft  [25,  21] 

■  LCV  ■ 
y  =  MCV 
NCV 

LCV  =  {l-U)LCVio  +  ULCVh, 

MCV  =  {l-fcv)MCV,o  +  UMCVM 
NCV  =  {l-fcv)NCV,o  +  fcvNCVH.i 
LCVhi  =  pcosa  +  rsina 
MCVhi  =  KqO  +  q 
NCVhi  =  KffP  —  /)sinQ  +  rcosa 
LCVio  =  p 
MCVio  =  q 
NCVio  =  r 

(3.16) 

where  fcv  is  a  blending  function  of  dynamic  pressure  that  blends  the  low  and  high  speed 
commands  as  shown  in  Fig. (3.3). 


Figure  3.3;  Command  variable  blending  function 


The  parameter  Ka  is  chosen  to  provide  a  natural  blend  of  pitch  and  vertical  acceleration 
since  the  longitudinal  stick  is  used  to  control  pitch  rotation  and  vertical  translation. 


Ka  = 


^s„,cz' 


mVc, 
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Table  3.1:  Command  variable  requirements  parameters 


Parameter 

Value 

^bJ 

qio 

19 

Qhi 

50^ 

a=0 


(3.17) 


Similarly,  the  parameter  Kp  'xs  chosen  to  provide  a  natural  blend  of  yaw  and  side  acceleration 
since  the  pedals  are  iised  to  control  yaw  rotation  and  side  translation. 


Ka  = 


iiS,e,c7j 


mVcc 


(3.18) 


Note  that  due  to  the  tailless  nature  of  the  ICE  aircraft,  CyJ  is  very  small  as  compared 

^  la=0 

to  similar  aircraft  with  vertical  tails  such  as  the  F-18  High  Angle-of-attack  Research  Vehicle 
(HARV) 


C 


O 


ICE 

a=0 
\HARV 

Q=0 


—.046  rad  ^ 
—1.2  rad“^ 


(3.19) 


UseofC^.;^=  J  in  eq.(3.18)  resulted  in  poor  sideslip  regulation.  This  led  to  the  most 

significant  tailless  aircraft  specific  dynamic  inversion  design  change  from  past  work  [25].  The 
choice  of  Kq  sho\ild  reflect  CyA  of  a  similar  vehicle  with  a  vertical  tail  otherwise,  poor 

^  ^)a=0 

performance  or  even  instability  may  result.  The  command  variable  requirements  parameters 
are  given  in  Table  3.1. 


Command  Shaping  Requirements 

The  choice  of  command  shaping  parameters  was  based  on  similar  aircraft  siich  as  the  F-16 
and  F-18.  The  initial  shaping  parameters  are  given  in  Table  3.2  and  will  be  tuned  during  the 
piloted  simulation  analysis  effort.  Note  that  variable  superscripts  refer  to  the  corresponding 
control  axis,  and  absence  of  a  superscript  implies  use  for  all  axes. 


Feedback  and  Feedforward  Augmentation  Requirements 

The  nominal  choice  of  aiigmentation  parameters  are  given  in  Table  3.3.  These  parameter 
choices  are  based  upon  uncertainty  bounds  and  flying  qualities  requirements  for  similar 
fighter  aircraft  [11,  25]. 
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Table  3.2:  Command  shaping  requirements  parameters 


Parameter 

JNominal  Value 

^roll 

.Ti 

ri,roll 

^roll 

Vi" 

rrToll 

vl°" 

rr'f'oll 

•''5 

!/r" 

pitch 

^^tch 

pitch 

X2 

^^tch 

vitch 
•'^3  , 

pitch 

2/3 

pitch 

X^ 

pitch 

Va 

vitch 
•'^5  ^ 

yaw 

x{ 

yaw 

y\ 

yaw 

X'2 

yaw 

y'2 

yaw 

•'^3 

yaw 

yl 

x^ 

yaw 

yl 

•^5 

yaw 

y'b 

LCVmax 

MCVmax 

NCVmax 

2.5  ^ 

sec 

0.0  lbs 

0.0  ^ 

sec 

1.0  lbs 

0.0  ^ 

sec 

4.0  lbs 

20.0  ^ 

sec 

8.0  lbs 

80.0  ^ 

sec 

16.57  lbs 

0.0  lbs 

0.0  ^ 

sec 

1.75  lbs 

0.0  ^ 

sec 

4.0  lbs 

5.0  ^ 

sec 

15.0  lbs 

40.0  ^ 

30.25  lbs 

0.0  lbs 

0.0  ^ 

sec 

5.0  lbs 

0.0  ^ 

sec 

96.0  lbs 

95.0  lbs 

95.0  lbs 
JVCV.^|g 
mm(|,-  200) 

100.0  ^ 

sec 

30.0  ^ 

sec 

Air  Force  Research  Laboratory 


20 


Table  3.3:  Augmentation  requirements  parameters 


Parameter 

JNominai  Value 

K  rad 
sec 

fi 

.25 

fc 

.5 

Table  3.4:  Allocation  requirements  parameters 


Parameter 

Description 

Nominal  Value 

AT 

Wd 

flight  control  system  update  rate 

control  deficient  weight 

0.01  sec 
r  1  0  0  1 
0  1  0 

0  0  1 

Control  Allocation  Requirements 

The  nominal  choice  of  control  allocation  requirements  parameters  are  given  in  Table  3.4. 

3.4.2  Configuration  Model 

The  configuration  model  contains  all  functions  and  parameters  that  pertain  to  the  aircraft 
configuration.  The  configuration  model  contains  an  airframe  dynamical  model,  actuator 
model  and  mission  segment  objective  models.  The  airframe  dynamical  model  provides  the 
stability  and  control  derivatives  to  the  augmentation  and  allocation  modules.  The  actuator 
model  provides  actuator  constraints  to  the  allocation  module.  The  mission  segment  objective 
models  provide  control  preferences  and  weights  for  control  sufficiency  optimization. 

The  airframe  model  required  the  most  development  effort,  and  it  still  needs  the  most 
refinement,  especially  at  high  speeds.  A  tradeoff  exists  between  airframe  model  fidelity  and 
controller  complexity.  A  high  fidelity  model  will  give  good  control  law  performance  but  will 
be  complex  requiring  significant  computational  throughput  and/or  memory.  On  the  other 
hand  a  low  fidelity  model  will  be  simple,  not  require  much  computational  throughput  and/or 
memory,  but  control  law  performance  may  be  poor.  Highly  robust  controllers  may  be  used  to 
account  for  low  model  fidelity,  but  may  be  so  complex  that  their  computational  throughput 
and  memory  requirements  may  negate  the  benefits  of  a  simple  model  [3]. 

Dynamic  inversion  is  chosen  as  the  flight  control  law  method  for  this  study  since  it  pro¬ 
vides  a  nice  compromise  between  controller  complexity  and  robustness  as  indicated  by  past 
studies  [22,  11,  25,  1,  21,  5].  The  airframe  model  is  developed  to  the  lowest  fidelity  possible 
such  that  closed-loop  stability  and  performance  is  adequate  as  determined  by  linear  analyt¬ 
ical  robustness  and  flying  qualities  tools  as  well  as  nonlinear  batch  and  piloted  simulations. 
This  lowest  fidelity  model  approach  is  taken  to  determine  the  level  of  model  fidelity  required 
by  the  dynamic  inversion  synthesis  method.  The  final  on-board  airframe  model  is  nonlinear, 
blit  a  scheduled  linear  model  and  a  perfect  model  were  developed  for  comparison  purposes. 
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Mach  number 


Figure  3.4:  Flight  Envelope  Linear  Model  Grid 

Scheduled  Linear  Airframe  Model 

A  full-envelope  airframe  model  consisting  of  linear  interpolation  between  linear  perturba- 
tional  models  was  initially  used  to  model  the  aircraft  dynamics.  Linear  models  were  gener¬ 
ated  at  the  flight  conditions  indicated  in  the  flight  envelope  in  Fig.(3.4). 

At  each  flight  condition,  the  elevens  and  pitch  flap  were  used  to  trim  the  aircraft  moments 
in  wings  level  flight  at  the  following  angles-of-attack 

o,,im  =  [  -2  0  4  8  12  16  20  24  28  f  deg.  (3.20) 

Central  difference  numerical  perturbations  of  the  nonlinear  aircraft  simulation  model  in 
eq.(2.1)  were  used  to  generate  10  state  linear  models  of  the  aircraft  rigid  body  dynamics  of 
the  following  form 

Xg  —  .AgXg  -h  BgU 

Xg  =  [uaq9h0pr(l>'tp]^ 

U  =  [  Sei  6pf  8(imtr  8pty  8ytv  ^o6//(  8oblfr  ] 

The  10  state  linear  models  were  reduced  to  8  state  linear  models  of  the  following  form  for 
control  design 

Xm  =  [aq^pr<j)eV]^ 
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Lateral  command  variable  (deg/sec) 


Figure  3.5:  Mach  0.35,  5000  ft.  Altitude.  Min  Vec  CA  Mode 


=  0 
=  0 
=  0 
=  0 
=  0 
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=  0 
=  0 

where  Fa,  Fg,  F^,  Fp,  Fj-  are  rows  of  F  in  eq.(2.1).  The  elements  of  the  and  matrices 
are  stored  in  data  tables  as  functions  of  Mach  nximber,  altitude,  and  angle-of-attack,  and 
linear  interpolation  is  used  between  data  points. 

As  the  simulation  analysis  proceeded,  it  became  clear  that  higher  fidelity  models  of  some 
effects  were  necessary.  The  command  variable  and  model  error  responses  in  Figs.(3.5)-(3.7) 
show  a  simulation  example  where  the  linear  model  is  inadequate.  The  command  variable 
tracking  with  the  linear  model  is  poor  relative  to  the  higher  fidelity  models  due  to  the  large 
model  errors. 
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Nonlinear  Airframe  Model 

A  nonlinear  on-board  model  was  designed  to  improve  control  law  performance  by  including 
nonlinear  state  dependencies  such  as  inertial  coupling  and  off-trim  effects  such  as  maneuver- 
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Directional  command  variable  (deg/sec) 


Figure  3.6:  Mach  0.35,  5000  ft.  Altitude,  Min  Vec  CA  Mode 


Longitudinal  control  error 


Figure  3.7:  Mach  0.35,  5000  ft.  Altitude,  Min  Vec  CA  Mode 
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ing  thrust.  Refinements  to  the  scheduled  linear  model  of  the  previous  section  were  made  to 
generate  the  nonlinear  model.  Least  squares  polynomial  data  fits  were  used  to  generate  some 
of  the  nonlinear  aerodynamic  and  propulsive  forces  and  moments.  The  nonlinear  on-board 
model  has  the  following  form 
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V,ssdr 

= 

R^^ 

V  ,ssdr 

jdNL 

^V.oblfl 

= 

jdLS 

^V,oblfl 

jdNL 

^V.oblfr 

jdLS 

^V.oblfr 

V 


u  — 

V  — 

w  = 

u  = 

y  =Z 

w  = 


y(r^M!an^ay'(l  +  tan^^ 

Vtan^ 


\/l 


+  tan^/3 
wtana 

qSrefCx  -\-T  —  DgCosa  —  H^sin^ 
r?.’  —  qw  - — - 


m 


pw  —  ru  + 


qSrefCy  +  VFcosfein^ 


m 


qu  —  pi>  + 


qSrefCz  +  Demina  +  W  cos0cos4> 


(3.22) 


The  aerodynamic  and  propulsive  parameters  from  eq.(3.22)  that  are  modeled  using  piece- 
wise  least  squares  polynomials  of  the  data  are  contained  in  Table  3.5.  Note  that  the  fit 
structure  in  Table  3.5  is  notional.  So  for  example,  although  the  coefficient  cq  appears  for 
many  parameters,  it  is  not  the  same  for  each  parameter. 
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Table  3.5:  Nonlinear  model  data  fits 


symbol 

T 


^f^-pflap 


c 

^Tlamt 


)escnption 

Thrust 


Engine  ram  drag 
X  body  axis  force  coefficient 
Y  body  axis  force  coefficient 
Z  body  axis  force  coefficient 


Roll  moment  coefficient 


Pitch  moment  coefficient 


Yaw  moment  coefficient 


Roll  moment  due  to  elevon 
Roll  moment  due  to 
all  moving  tip 
Roll  moment  due  to 
spoiler  slot  deflector 
Roll  moment  due  to 
outboard  leading  edge  flap 
Pitch  moment  diie  to 
elevon 

Pitch  moment  due  to 
pitch  flap 

Pitch  moment  due  to 
all  moving  tip 
Pitch  moment  due  to 
spoiler  slot  deflector 
Yaw  moment  due  to 
all  moving  tip 
Yaw  moment  due  to 
spoiler  slot  deflector 
Yaw  moment  due  to 
outboard  leading  edge  flap 


Pit  Structure  _ 

Co  +  CiM  +  C2M^ 

+czM^  +  h  (c4  +  C5M 

+C6A/^  +  cjM^) 

Co  +  CiM  +  C2M^ 

+/?.(C3  +  C4A'/  + 

Co  +  cia  +  C2Q;^ 

+P‘^{c3  +  QCk  +  +  cea^) 

0{co  +  CiQ  +  C2a^  +  +  C4a:‘‘) 

+  ^(C5P+C6r) 

Co  +  CiQ  +  020^  +  Csa^ 

+/?^(C4  +  C5Q  +  Cea^  +  c^a^)  +  cs^ 

/?(co  +  cia  +  C2Q^  +  +  C4a^  -f  c^a^) 

+  w  ((*^0  +  +  ‘^3^^  +  C4tt'‘)p 

+  (co  +  CjO  +  C20'^  +  +  C4Q:^)r) 

(co(a)  +  ci(Q:)a  +  C2(a)a;^ 

+C3 (a)Q^)  +  (co  +  Cia  +  C2a^  _ 
+C3Q^  +  C4a'‘  +  C5M  +  cqMo)  ^ 

Co  +  CiO  +  C2a;^  +  Czo? 

+/?^(c4  +  c-^Oi  +  ceo;^  +  cja^) 

+W  ((^0  +  CiQ  +  C2a^  +  C3a^  +  C4a'')7; 
+(co  +  cja  +  C2a^  +  C3a^  +  C4a^)r) 

Co  +  CiQ  +  020^  +  csa^ 

Co  +  CiO  +  C2a^  +  C3a^ 

Co  +  ciO  +  C2a^  +  cstt^ 

Co  +  CiQ  +  C2a^  +  CsO:^ 

Co  +  CiQ  +  C2a^  +  C3Q;^ 

Co  +  Cia  +  C2a^  +  C3a^ 

co(a)  +  Cl  (a)a  +  C2(a)a^  +  C3{a)a^ 


Co  +  CiO  +  C2a^  +  csa^ 
Co  +  ciQ  +  C2a^  +  c^a^ 
Co  +  C'lQ  +  020^  +  C3Q;^ 
Co  +  Cio  +  C2a^  +  Cza^ 
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Figure  3.8:  Static  pitching  moment  spline  fit 


Although  not  explicitly  annotated,  the  propulsion  parameter  (T,  Dg)  coefficients  are 
functions  of  power  lever  angle.  Recall  that  the  simulation  aerodynamic  model  has  a  low 
speed  (Mach<  0.5)  and  a  high  speed  part  (Mach>  0.6)  and  are  blended  together  at  inter¬ 
mediate  speeds  (0.5<Mach<0.6).  The  aerodynamic  parameter  coefficients  in  Table  3.5  are 
also  functions  of  these  Mach  partitions. 

The  dynamic  inversion  controller  was  particularly  sensitive  to  the  low  speed  aerodynamic 
static  pitching  moment,  so  high  fidelity  cubic  splines  were  used  as  a  model.  A  polynomial 
representation  of  the  static  pitching  moment  cubic  spline  is  given  from  Table  3.5  as 

Cmo  =  co(a)  -I-  Cl (a)a;  +  C2(a)Q^  +  €3(0)0^  (3.23) 

The  final  spline  fit  is  shown  in  Fig.(3.8). 


Perfect  Airframe  Model 

Besides  a  linear  and  nonlinear  on-board  aircraft  model,  a  “perfect”  on-board  model  mode 
was  implemented  for  comparison  purposes.  This  “perfect”  on-board  model  is  truly  only 
nearly  perfect  since  it  has  the  following  affine  form 


Xr, 


=  A 


PERF 


+  B.^ 


PERF 


U 
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aPERF 


pL 

Op 

Of 

0 


(3.24) 


where  aa,Oq,ap,af  are  elements  of  a  in  eq.(2.2)  and  is  identical  to  in  eq.(3.22) 

except  that  the  true  thrust  T  and  ram  drag  data  are  used.  Of  course  this  mode  is  not 
intended  to  be  a  viable  alternative  for  implementation,  but  it  does  provide  a  useful  mode  for 
comparison. 


Actuator  model 

The  on-board  actuator  model  consists  of  actuator  position  and  no-load  rate  limits.  The 
hinge  moment  effect  for  all  surfaces  were  neglected  since  the  simulation  model  also  did  not 
include  these  effects.  Actuator  dynamics  are  neglected  since  the  control  allocation  approach 
is  inherently  static.  Control  effector  rate  and  position  limits  are  derived  from  the  simulation 
model  data  in  Table  2.1  and  are  given  by  the  following 

=  [30  30  30  60  60  10  10  10  10  40  40  Y 
=  [  -30  -30  -30  0  0  -10  -10  0  0  0  0  Y 
Ur  =  [  150  150  50  150  150  60  60  150  150  40  40  Y  (3.25) 

Note  that  for  the  on-board  model,  the  spoiler  slot  deflector  (elements  8,9)  upper  limits  are 
artifically  redticed  from  the  actual  actuator  data  in  Table  2.1.  This  was  done  to  limit  the 
spoiler  slot  deflector  and  elevon  interactions. 

Mission  Segment  Objective  Model 

The  following  models  are  used  by  the  control  allocation  algorithm  to  optimize  specific  mission 
segment  objectives  such  as  minimum  drag  or  wing  loads.  Recall  that  a  mission  segment 
objective  is  completely  specified  by  a  control  preference  (wpre/)  corresponding  weight 
(VTu)  from  eq.(3.12). 

Minimum  Control  Deflection  The  minimum  control  deflection  control  allocation  mode 
is  implemented  by  the  following  choice  of  parameters 

Upref  =  [00000000000  ]^ 

=  [  1  1  1  1  1  1  1  1  1  1  1  ]  (3.26) 

which  equally  emphasizes  all  effectors  toward  zero  deflections. 
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Figure  3.9:  Elevon  minimum  drag  deflection 

Minimum  Drag  At  low  subsonic  speeds,  Mach<  0.5,  the  minimum  drag  control  allocation 
mode  uses  an  on-board  drag  model  to  compute  Upref-  Due  to  time  limitations,  no  on-board 
drag  model  was  developed  above  Mach  0.5.  So  for  Mach>  0.5,  the  minimum  drag  control 
allocation  mode  is  exactly  the  same  as  the  minimum  deflection  mode. 

The  low  speed  on-board  drag  model  was  generated  from  the  control  surface  axial  force 
wind  tunnel  data.  The  minimal  axial  force  deflections  for  the  elevons,  pitch  flaps,  and 
outboard  leading  edge  flaps  are  given  in  Figs. (3.9,  3.10,  3.11)  respectively.  The  minimum 
drag  deflection  data  is  questionable  since  large  constant  deflections  result  in  the  minimum 
axial  forces  which  seems  counter  intuitive.  Although  the  wind  tunnel  data  looks  suspect, 
the  minimum  drag  control  allocation  mode  may  be  demonstrated  using  this  data.  As  better 
data  is  collected,  the  on-board  model  may  also  be  revised. 

The  minimum  drag  control  allocation  mode  is  implemented  by  the  following  choice  of 
parameters  for  Mach  <  0.5 

Upref  =  [  felia)  fel{oi)  fpf{Oi)  0  0  0  0  0  0  foblf{<^)  foblf{o^)  f 

Wu  =  [  10.0  10.0  10.0  0.1  0.1  0.01  0.01  100.0  100.0  1.0  1.0  ]  (3.27) 

where  fei{oi),  fpf{a),  fobif{ot)  are  the  on-board  drag  due  to  control  deflection  models  indicated 
in  Figs.(3.9,  3.10,  3.11).  The  minimum  drag  control  allocation  mode  is  implemented  by  the 
following  choice  of  parameters  for  Mach  >  0.5 

=  [  0  0  0  0  0  0  0  0  0  0  0 

V14  =  [  1  1  1  1  1  1  1  1  1  1  1  ]  (3.28) 

which  are  identical  to  the  minimum  control  deflection  values. 
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Figiire  3.10:  Pitch  flap  minimum  drag  deflection 


Figure  3.11:  Outboard  leading  edge  flap  minimum  drag  deflection 
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Minimum  Wing  Loading  The  minimum  wing  load  model  is  simplistic.  Higher  weights 
are  chosen  for  outboard  control  effectors  since  they  qualitatively  give  higher  wing  root  bend¬ 
ing  moments  than  inboard  surfaces.  The  minimum  wing  loading  control  allocation  mode  is 
implemented  by  the  following  choice  of  parameters 

Mpre/  =  [00000000000]^ 

Wu  =  [  1.0  1.0  0.1  1000.  1000.  0.01  0.01  1.0  1.0  10.0  10.0  ]  (3.29) 

which  emphasize  the  all  moving  tips  toward  a  zero  deflection. 


Minimum  Signature  The  minimum  radar  signature  model  is  also  simplistic.  Higher 
weights  are  chosen  for  spoiler  slot-deflectors  and  all-moving  tips  since  they  qualitatively  give 
higher  radar  signatures  than  the  other  effectors.  The  minimum  signature  control  allocation 
mode  is  implemented  by  the  following  choice  of  parameters 

Uj,ref  =  [00000000000]^ 

W„  =  [  1.0  1.0  1.0  10.  10.  0.01  0.01  100.0  100.0  1.0  1.0  ]  (3.30) 


Null  Space  Injection  The  only  parameter  that  is  required  for  the  null  space  injection 
mode  is  the  nominal  weight  W  in  eq.(3.14)  chosen  to  equally  distribute  commands.  The 
elements  of  W  were  nominally  chosen  as  the  inverse  of  the  absolute  value  of  the  corresponding 
effector  deflection  limit.  These  were  tuned  to  the  following  values  based  on  simulation 
responses  to  achieve  equally  distributed  commands. 


\30’  30’  30’  120’  120’  15’  15’  120’  120’  80’  80J 


(3,31) 


Minimum  Thrust  Vectoring  The  minimum  thrust  vectoring  control  allocation  mode  is 
implemented  by  the  following  choice  of  parameters 

uj^ef  =  [00000000000]^ 

=  [  1.0  1.0  1.0  1.0  1.0  1000.0  1000.0  1.0  1.0  1.0  1.0  ]  (3.32) 

which  emphasize  the  thrust  vector  nozzle  toward  zero  deflection. 


Constants  and  Measurements 

Required  constants  for  the  configuration  model  are  contained  in  Table  3.6,  and  required 
measured  variables  are  contained  in  Table  3.7. 
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Table  3.6:  Nonlinear  model  constants 


Symbol 

Description 

Value 

^La  L  =  o 

lift  coefficient  due  to  angle  of  attack 
at  zero  angle  of  attack 

2-42;^ 

Ko 

crossover  velocity 

Ixx 

■T-axis  moment  of  inertia 

35, 479slug  -  ff 

^yy 

y-axis  moment  of  inertia 

78,  ATilslug  —  ff^ 

hz 

2-axis  moment  of  inertia 

110, 627slug  —  ft^ 

Ixz 

.T2-axes  product  of  inertia 

-525slug  -  ft^ 

Sref 

wing  planform  area 

808.6/^2 

c 

wing  mean  aerodynamic  chord 

28.75/t 

b 

wing  span 

37.5/t 

m 

aircraft  mass 

1017.9  slug 

W 

aircraft  gross  weight 

32, 750/5 

drwL 

waterline  distance  between  center  of  gravity 
and  thrust  application  point 

0.502b  ft 

drFS 

fuselage  station  distance  between  center  of  gravity 
and  thrust  application  point 

18.76// 

dowL 

waterline  distance  between  center  of  gravity 
and  ram  drag  application  point 

-0.3308// 

doFS 

fuselage  station  distance  between  center  of  gravity 
and  ram  drag  application  point 

-12.66// 

Table  3.7;  Non: 


inear  model  required  measurements 


Symbol 

Description 

Units 

V 

roll  rate 

deg 

sec 

Q 

pitch  rate 

deg 

sec 

r 

yaw  rate 

deg 

a 

angle  of  attack 

deg 

0 

angle  of  slip 

deg 

V 

total  velocity 

IL 

sec 

Q 

dynamic  pressure 

lb 

0 

roll  attitude  angle 

deg 

e 

pitch  attitude  angle 

deg 
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Summary 

The  on-board  model  contains  configuration  and  requirements  data.  Configuration  data  in¬ 
cludes  aerodynamic,  propulsion  and  actuator  models.  Requirements  data  includes  command 
variable  definitions  and  desired  command  variable  band  widths.  The  configuration  data  has 
the  following  known  modeling  errors  compared  to  the  simulation  model:  neglected  actuator 
dynamics,  aerodynamic  and  propulsive  approximations,  control  effector  interaction  effects 
neglected,  velocity  changes  neglected. 


3.5  Model  Error  Compensation 

The  on-board  model  in  Section  3.4  will  contain  errors  from  many  possible  sources.  The 
robustness  of  the  control  law  accounts  for  some  level  of  model  error.  Adaptation  may  be 
used  to  account  for  higher  levels  of  model  error  due  to  failures,  damage  or  highly  nonlinear 
phenomena  [7,  12,  19].  Indirect  and  direct  adaptation  methods  have  already  been  success¬ 
fully  demonstrated  on  flight  test  aircraft  [2,  5].  Future  plans  include  integration  of  on-line 
system  identification  [13]  and  adaptive  neural  networks  [28]  to  provide  adaptive  model  error 
compensation. 
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Chapter  4 

Flight  Control  Law  Analysis 


Analysis  of  the  flight  control  law  in  Chapter  3  is  contained  in  this  chapter.  Eigenvalue  anal¬ 
ysis  is  used  to  analyze  the  nominal  stability,  flying  qualities  criteria  are  used  for  nominal 
performance  analysis  and  structured  singular  values  and  batch  simulations  are  used  to  an¬ 
alyze  robust  stability  and  performance.  Real  time  simulation  analysis  has  begun  to  assess 
pilot  in  the  loop  performance,  and  results  will  be  reported  in  a  future  document. 

Eigenvalue  and  batch  simulation  analyses  are  presented  at  four  representative  flight  con¬ 
ditions.  The  flrst  flight  condition  is  Mach  0.35  and  5,000  ft.  altitude  and  represents  a  low 
dynamic  pressure  powered  approach  condition.  The  second  flight  condition  is  Mach  0.60  and 
15,000  ft.  altitude  and  represents  a  center  of  the  envelope,  air  combat  maneuvering  condi¬ 
tion.  The  third  flight  condition  is  Mach  0.90  and  1,000  ft.  altitude  and  represents  a  low 
altitude  high  subsonic  ingress  condition.  The  final  flight  condition  is  Mach  1.20  and  22,500 
ft.  altitude  and  represents  a  supersonic  cruise  condition.  Table  4.1  lists  the  properties  of 
these  four  analysis  flight  conditions. 

Structured  singular  value  and  flying  qualities  analyses  are  presented  for  a  subset  of  the 
flight  conditions  in  Fig. (3.4).  These  analysis  flight  conditions  provide  a  complete  dynamic 
pressure  variation  of  the  subsonic  flight  envelope  and  are  given  in  Table  4.2. 


Table  4.1:  Eigenvalue  and  batch  simulation  analyses  flight  cond 


Description 

Motation 

Mach 

Altitude 

(ft) 

Angle  oi 
attack  (deg) 

Dynamic 
pressure  (^) 

Low  q,  powered  approach 

m35h5 

0.35 

8.01 

151 

Envelope  center,  air  combat 

m6hl5 

■inuB 

3.70 

High  subsonic,  ingress 

m9hl 

HnMlyyH 

0.904 

1157 

Supersonic,  cruise 

■Jail 

1.48 

882 

itions 
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4.1  Nominal  Stability 


Nominal  stability  of  the  closed  loop  flight  control  system  is  guaranteed  from  the  dynamic 
inversion  flight  control  theory  if  the  chosen  command  variable  dynamics  are  stabilizing  and 
the  aircraft  is  minimum  phase  (stable  zero  dynamics),  see  Section  2.2.  The  stability  of  the 
command  variable  and  zero  dynamics  is  determined  by  eigenvalue  analysis  and  confirmed  by 
simulations. 

The  nominal  closed  loop  command  variable  dynamics  are  linear  by  design  and  given  by 


y  = 

T-i 


0  u;c  0  0 

0  0  0 
—Uc  0  0  tJc 

0  0  0  0 

0  0  0  0 

-ojji  0  0  0 


of  the  following  system  matrix 

-UJc  0 

0  —^c 

,  _  0  0 
^cv  =  0 

0  -Uc] 

0  0 


eig(A„;)  = 


—Uc  0  0  Wc  0  0 

0  —Uc  0  0  LJc  0 

0  0  — a;c  0  0 

-oJcfi  0  0  0  0  0 

0  -Ucfi  0  0  0  0 

0  0  -Ucfi  0  0  0 

-2.5,  -2.5,  -2.5,  -2.5,  -2.5,  -2.5 


^cfc 

0 

0  1 

0 

i^cfc 

0 

0 

0 

(^cfc 

0 

0 

0 

^cfi 

0 

0 

0 

- 

?ly  specified  by  the 

-5 

0 

0  5 

0 

-5 

0  0 

0 

0 

-5  0 

-5i 

4 

0 

0  0 

0 

0  0 

0 

0 

1 

o 

which  proves  nominal  asymptotic  stability  of  the  command  variables  since  all  eigenvalues 
have  negative  real  parts. 

The  stability  analysis  of  the  zero  dynamics  is  more  complicated  due  to  nonlinearities. 
The  assessment  of  nominal  stability  is  made  by  analyzing  the  following  linear  approximation 
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to  the  complementary  dynamics  in  eq.(2.8) 

a(2)  =  -W-'B^  {B,W-'Biy' A,,z  (4.3) 

which  assumes  the  following  linear  form  for  the  control  allocator 

tx  =  -W-^Bl  (ByW-^Biy^  (4.4) 

The  asymptotic  stability  of  the  system  in  eq.(4.3)  is  determined  by  the  eigenvalues  of  the 
following  system  matrix 

Azd  =  A,,  -  B,W-^Bl  {ByW-^B^y'  Ay,  (4.5) 

The  linear  control  allocation  weight  in  eqs.(4.3,  4.4,  4.5)  is  chosen  to  best  reflect  simulation 
responses  of  the  linear  programming  control  allocation  method  in  Section  (3.3).  Simulation 
responses  have  shown  that  the  linear  programming  control  allocation  method  attempts  to  put 
as  many  control  effectors  on  their  preferences  as  possible  thus  only  using  a  few  to  achieve 
the  control  law  moment  command  automatically  canceling  the  residual  moment  effect  of 
those  effectors  on  their  preference.  Recall  that  most  control  allocation  modes  in  Section 
(3.3)  have  zero  preferences.  Thus,  the  linear  control  allocation  function  in  eq.(4.4)  is  a  good 
approximation  to  the  linear  program  control  allocator  if  very  large  weights  corresponding 
to  those  effectors  that  follow  zero  are  chosen.  Weights  are  chosen  to  approximate  the  linear 
program  control  allocation  at  the  corresponding  flight  condition  based  on  simulation  time 
histories  of  the  control  deflections.  For  example,  the  linear  program  primarily  uses  the 
elevens  and  all  moving  tips  for  small  inputs  at  Mach=0.35  and  5000  ft.  altitude.  Therefore, 
all  diagonal  elements  of  the  inverse  weight  matrix  at  this  flight  condition  are  chosen  zero 
except  those  corresponding  to  the  elevens  (elements  1  and  2)  and  all  moving  tips  (elements 
4  and  5).  Similar  weights  are  constructed  based  on  simulation  responses  for  the  other  flight 
conditions  as  well.  Two  other  weights  were  generated  for  illustrative  comparison  purposes. 
One  weight  was  constructed  to  mimic  a  control  allocator  that  only  uses  spoiler  slot  deflectors 
for  yaw  control.  The  final  weight  was  an  identity  matrix  which  corresponds  to  tising  all 
surfaces.  The  fomr  weights  used  for  analysis  are  given  by  the  following 

=  diag(|l  1  0  1  1  0  0  0  0  0  0  ]) 

=  diagdl  1  1  0  0  1  1  0  0  0  0]) 

=  <iiag([l  1  0  0  0  1  1  0  0  0  0]) 

M'm.W  =  <liag([l  1  0  0  0  0  1  0  0  0  0]) 

W'r.J  =  diagdl  1  1  0  0  1  0  1  1  0  0]) 

Wf^,  =  diag([  11111111111])  (4.6) 

Table  4.3  shows  the  eigenvalues  of  the  zero  dynamics  system  matrix  in  eq.(4.5)  at  each 
analysis  flight  condition  in  Table  4.1  for  each  weight  in  eq.(4.6).  Note  that  the  diagonal 
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Table  4.3:  Subroutine  Functionality 


mi 

■smf 

mmumamm 

MilMIBIlliKflM 

-7.46,  -2.72 

-5.24,  -2.86 

BIEWBilliWW 

■WEiaroiBM 

mmmm 

mmammsftm 

12/7.225 

mMYmmmm 

mmmmsim 

■rgBMaiBM 

^  ssd 

mmxmtwifam 

^  ident 

mMYmamm 

mrKiimmum 

elements  of  Table  4.3  correspond  to  the  linear  approximation  of  the  nominal  allocation  mode 
for  the  particular  flight  condition. 

It  is  seen  that  the  stability  of  the  zero  dynamics  is  relatively  insensitive  to  the  choice  of 
weights  at  all  flight  conditions  except  the  low  dynamic  pressure  condition.  Although  Wrn^bhb 
gives  nominally  stable  zero  dynamics  at  Mach  0.35  and  5000  ft.  altitude,  W^sd  destabilizes 
the  zero  dynamics.  Small  roll  stick  pulse  command  simulation  responses  in  Figs.  (4. 1,4.2) 
confirm  the  eigenvalue  stability  analysis  of  the  zero  dynamics  at  this  flight  condition.  It  is 
seen  that  the  aircraft  departs  when  using  only  spoiler  slot  deflectors.  Note  that  SSD  refers 
to  the  control  allocation  that  primarily  uses  the  spoiler  slot  deflectors  for  yaw  control,  and 
AMT  refers  to  the  control  allocation  that  primarily  uses  the  all  moving  tips  for  yaw  control. 
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Lateral  command  variable  (deg/sec) 


Figure  4.1:  Mach  0.35,  5000  ft.  Altitude,  Zero  Dynamics  Stability 


This  analysis  shows  that  the  linear  program  control  allocation  method  may  reduce  sta¬ 
bility  margins  or  even  destabilize  the  zero  dynamics  when  it  does  not  use  all  of  the  effectors. 
This  possibly  explains  the  lightly  damped  behavior  in  the  simulation  responses  for  some 
control  allocation  modes  at  Mach  0.35  and  5000  ft.  altitude  in  Section  (4.4). 


4.2  Flying  qualities 

The  MATLAB  flying  qualities  toolbox  [14]  was  used  to  initially  analyze  the  flying  qualities  of 
the  flight  control  system  at  the  flight  conditions  in  Table  4.2.  The  longitudinal  flying  qualities 
were  analyzed  using  the  bandwidth  criterion  and  the  control  anticipation  parameter  (CAP) 
from  a  low  order  equivalent  system  (LOES).  It  is  seen  in  Fig.  (4.3)  that  Level  1  CAP  is 
achieved  and  borderline  Level  1,  Level  2  bandwidth  is  achieved.  The  lateral/directional 
flying  qualities  were  analyzed  using  LOES  parameters  and  the  roll  pilot  induced  oscillation 
(PIO)  criterion.  It  is  seen  in  Figs.(4.4,4.5)  that  Level  1  flying  qualities  are  achieved  and  PIO 
tendencies  do  not  exist. 
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Figure  4.4:  Lateral/Directional  Flying  Qualities  Criteria 


More  details  of  this  analysis,  as  well  as  a  comparison  with  a  linear  parameter  varying 
control  law,  appear  in  [9]. 


4.3  Robustness 


Structured  singular  values  were  used  to  assess  robust  stability  and  performance  of  the 
flight  control  system  at  the  flight  conditions  in  Table  4.2.  Separate  longitudinal  and  lat¬ 
eral/directional  analyses  were  performed.  Note  that  identity  weights  were  used  for  the  linear 
control  allocator  approximation  in  eq.(4.4)  for  this  analysis.  Figure  (4.6)  shows  the  struc¬ 
tured  singular  value  analysis  model  for  the  lateral/directional  axes.  There  are  three  un¬ 
certainty  blocks.  The  first  uncertainty  block  is  unstructured  with  2  inputs  and  2  outputs. 
Command  response  performance  is  captured  in  this  block  with  the  following  roll  and  yaw 
error  weight 

w  = - 1 - 

®  .01s2  + 0.14s +  1 

which  is  a  low  pass  filter  to  emphasize  low  frequency  tracking.  The  second  uncertainty  block 
is  unstructured  with  2  inpiits  and  2  outputs.  Controller  bandwidth  constraints  are  captured 
in  this  block  with  the  following  roll  and  yaw  input  uncertainty  weight 


Wi  = 


10s -H 
s  +  100 


(4.8) 


which  a  high  pass  filter  with  a  gain  crossover  frequency  of  10^  to  emphasize  high  frequency 
attenuation.  The  third  uncertainty  block  is  a  structured  uncertainty  block  with  2  inputs 
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Figure  4.5:  Lateral/Directional  Flying  Qualities  Criteria 


and  2  outputs.  Stability  derivative  parameter  uncertainty  of  0.2N^  and  0.5L^  is  captured  in 
this  block.  The  lateral/directional  structured  singular  values  for  robust  stability  and  robust 
performance  are  shown  in  Fig. (4.7).  Since  the  maximum  structured  singular  value  for  robust 
stability  is  approximately  1.4  for  all  flight  conditions,  robust  stability  is  guaranteed  at  all 
flight  conditions  for  r^Aiat/dir  where  Aiat/dir  represents  the  uncertainty  in  Fig.(4.6).  Similarly 
since  the  maximum  structured  singular  value  for  robust  performance  is  approximately  1.9  for 
all  flight  conditions,  robust  performance  is  guaranteed  at  all  flight  conditions  for  -^Aiat/dir- 
Figure  (4.8)  shows  the  structured  singular  value  analysis  model  for  the  longitudinal  axis. 
There  are  three  uncertainty  blocks.  The  first  uncertainty  block  is  unstructured  with  1  input 
and  1  output.  Command  response  performance  is  captured  in  this  block  with  the  error 
weight  in  eq.(4.7).  The  second  uncertainty  block  is  unstructured  with  1  input  and  1  output. 
Controller  bandwidth  constraints  are  captured  in  this  block  with  the  input  uncertainty  weight 
in  eq(4.8).  The  third  uncertainty  block  is  a  structured  uncertainty  block  with  2  inputs 
and  2  outputs.  Stability  derivative  parameter  uncertainty  of  0.1  Zq  and  I.OMq  is  captured 
in  this  block.  The  longitudinal  structured  singular  values  for  robust  stability  and  robust 
performance  are  shown  in  Fig. (4.9).  Since  the  maximum  structured  singular  value  for  robust 
stability  is  approximately  0.9  for  all  flight  conditions,  robiist  stability  is  guaranteed  at  all 
flight  conditions  for  ^ Along  where  Along  represents  the  uncertainty  in  Fig.  (4.8).  Similarly 
since  the  maximum  structured  singular  value  for  robust  performance  is  approximately  1.2 
for  all  flight  conditions,  robust  performance  is  guaranteed  at  all  flight  conditions  for  -^Aiong- 
More  details  of  this  analysis,  as  well  as  a  comparison  with  a  linear  parameter  varying  control 
law,  appear  in  [9]. 
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Figure  4.7:  Lateral/Directional  Axes  Structured  Singular  Values 


4.4  Batch  simulation 

The  non-real  time  batch  simiilation  analysis  was  done  using  the  GENESIS  simulation  en¬ 
vironment  to  assess  robust  stability  and  performance  at  the  flight  conditions  in  Table  4.1. 
There  are  two  parts  to  the  batch  simulation  analysis:  control  allocation  mode  verification 
analysis  and  final  control  law  analysis.  The  objectives  of  the  control  allocation  mode  verifica¬ 
tion  analysis  were  to  verify  that  each  mode  minimized  the  corresponding  objective,  identify 
any  problems  associated  with  each  mode  and  choose  nominal  control  allocation  modes  for 
each  flight  condition.  The  objective  of  the  final  control  law  analysis  was  to  evaluate  the 
entire  control  law  with  the  nominal  control  allocation  mode  choice  at  each  flight  condition. 

The  presentation  of  simulation  time  histories  is  consistent  throughout  this  document. 
The  inner  loop  command  variable  responses  are  plotted  with  the  corresponding  reference 
commands  (CMD).  Left  and  right  control  surface  deflections  (DEF)  are  plotted  with  the 
corresponding  actuator  command  (CMD)  and  control  allocation  optimal  preference  (PREF). 
Symmetric  and  asymmetric  deflections  are  also  plotted  to  show  contributions  of  control  pairs 
to  each  axis.  Control  law  command  errors  are  plotted  which  indicate  when  there  is  insuf¬ 
ficient  control  power  to  achieve  the  control  law  command.  Errors  in  the  on-board  models 
of  the  aerodynamics  and  actiiators  are  indicated  for  each  axis.  Aerodynamic,  specifically 
control  derivative,  model  errors  are  indicated  by  differences  between  model  and  truth  accel¬ 
eration  responses.  Actuator  model  errors,  specifically  due  to  neglecting  actuator  bandwidth 
limitations,  are  indicated  by  differences  between  control  law  commands  (CMD)  and  achieved 
accelerations  due  to  control  deflections  (OUT). 
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Figure  4.9:  Longitudinal  Axis  Structured  Singular  Values 


4.4.1  Control  Allocation  Mode  Verification  Analysis 

The  control  allocation  mode  verification  analysis  was  done  for  a  single  loaded  roll  maneuver 
at  a  single  flight  condition  for  all  control  allocation  modes.  The  low  dynamic  pressure  flight 
condition  at  Mach  0.35  and  5000  ft.  altitude  was  chosen  for  this  analysis  since  the  zero 
dynamics  at  this  flight  condition  were  the  least  stable  and  were  most  sensitive  to  control 
allocation  modes  as  indicated  in  Sec.  4.1. 

Performance  of  the  control  allocation  modes  is  evaluated  by  metrics  that  correspond  to 
the  intended  objective  of  each  mode.  The  following  control  allocation  metrics  were  used  for 
evaluation 
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Figure  4.10:  Control  Allocation  Mode  Metrics 
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where  AT  is  the  simulation  time  step,  t/  is  the  simulation  run  time,  6{k)  is  the  vector  of 
control  deflections  at  the  simulation  time  step,  and  Coik)  is  the  coefficient  of  drag.  The 
regressor  matrix  H{k)  consists  of  symmetric  left  and  right  control  pairs  and  pitch  thrust 
vectoring  for  the  longitudinal  axis  analysis.  It  consists  of  asymmetric  left  and  right  control 
pairs  and  yaw  thrust  vectoring  for  the  lateral/directional  axes  analysis.  The  H{k)  matrix 
has  data  from  the  previous  0.8  sec  or  80AT,  and  each  column  is  scaled  by  the  maximum 
absolute  value  of  the  column. 

The  functionality  of  the  control  allocation  modes  is  verified  by  the  bar  charts  in  Figs.  (4.10)- 
(4.12).  Each  figure  corresponds  to  an  evaluation  metric  in  eq.(4.9),  and  each  bar  in  the  figures 
corresponds  to  a  control  allocation  mode.  For  example,  the  second  bar  in  the  first  figure  is  the 
minimum  drag  control  allocation  mode  performance  measured  by  the  minimum  deflection 
metric  Jdej- 

It  is  seen  that  all  control  allocation  modes  minimize  the  corresponding  objective  except 
minimum  drag  mode.  All  of  the  modes  give  approximately  the  same  drag.  This  is  due  to  the 
fact  that  drag  due  to  control  deflections  is  a  very  small  part  of  total  drag.  Therefore,  any 
drag  benefits  from  minimum  drag  control  allocation  does  not  seem  to  be  worth  the  added 
complexity  of  modeling  control  effector  drag  in  the  control  law.  It  is  interesting  to  note  that 
the  null  space  injection  mode  increases  the  deflection  metric  by  only  56%  while  minimum  drag 
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Figure  4.11:  Control  Allocation  Mode  Metrics 


Thrust  Vector  metric 
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mode  increases  the  deflection  metric  by  503%.  This  shows  that  extremely  large  deflections  do 
not  always  increase  drag.  The  minimum  wing  load,  minimum  signature  and  minimum  vector 
modes  resist  motion  of  specific  surfaces  which  result  in  very  poor  performance  of  modes  that 
use  those  surfaces.  The  minimum  deflection  mode  does  well  for  each  metric  as  it  is  not  the 
worst  performing  mode  for  any  metric.  Another  interesting  note  is  that  the  injection  mode 
performs  similarly  to  the  minimum  deflection  mode  except  for  minimizing  signature. 

The  purpose  of  the  injection  mode  is  to  enhance  on-board  system  identification  for  control 
law  reconfiguration  capability.  Although  no  on-line  system  identification  was  included  in  the 
flight  control  system  of  the  present  study,  it  could  easily  be  integrated  into  the  modular 
control  law  architecture  to  enhance  the  fidelity  of  the  parameters  from  the  on-board  model. 
It  has  been  shown  that  the  condition  number  of  a  matrix  containing  scaled  control  deflections 
is  a  good  indicator  of  system  identification  performance  [29].  Low  condition  numbers  may 
result  in  good  parameter  estimates,  while  high  condition  numbers  will  certainly  result  in 
poor  estimates.  It  has  also  been  shown  [29]  that  condition  numbers  less  than  250  will 
typically  give  good  estimates.  The  inverse  of  the  condition  numbers  of  a  matrix  containing 
scaled  control  deflections  are  plotted  in  Figs. (4. 13)- (4. 18).  Note  that  larger  numbers  indicate 
better  conditioned  matrices  since  condition  number  inverses  are  plotted.  The  straight  line  in 
all  the  plots  correspond  to  a  condition  number  of  250.  It  is  seen  that  all  the  control  modes 
result  in  good  condition  numbers  only  during  maneuvering  except  the  injection  mode  which 
results  in  good  condition  numbers  all  the  time. 
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Figure  4.13:  Mach  0.35,  5000  ft.  Altitude,  Min  Def  CA  Mode 
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Injection  metric  —  lat/dir 


Figure  4.14:  Mach  0.35,  5000  ft.  Altitude,  Min  Drag  CA  Mode 


Figure  4.15:  Mach  0.35,  5000  ft.  Altitude,  Min  Load  CA  Mode 
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Figure  4.16:  Mach  0.35,  5000  ft.  Altitude,  Min  Signature  CA  Mode 


Figure  4.17:  Mach  0.35,  5000  ft.  Altitude,  Injection  CA  Mode 
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Injection  metric  -  lat/dir 


Figure  4.18;  Mach  0.35,  5000  ft.  Altitude,  Min  Vec  CA  Mode 


The  plots  in  Figs.(4.19)-(4.84)  are  batch  simulation  time  histories  of  a  loaded  roll  ma¬ 
neuver  at  the  low  dynamic  pressure  condition  for  the  5  control  allocation  modes. 

Minimum  deflection  mode 

The  minimum  deflection  control  allocation  mode  time  histories  are  contained  in  Figs.  (4.19)- 
(4.29).  Note  that  there  is  only  slight  coupling  between  command  variable  axes,  and  angle  of 
slip  regulation  is  very  good.  Elevons  and  all  moving  tips  are  the  most  active  control  effectors 
with  all  other  surfaces  only  deviating  from  their  preferences  at  the  onset  of  stick  commands 
when  insufficient  control  power  exists  to  achieve  the  control  law  acceleration  commands. 
Note  the  slight  directional  oscillations  prevalent  in  the  yaw  thrust  vector  and  directional 
command  variable  responses.  This  is  due  to  the  lower  bandwidth  actuator  of  the  thrust 
vector  nozzle  compared  to  the  aerodynamic  effectors.  There  is  almost  no  on-board  actuator 
model  errors  and  only  slight  aerodynamic  model  errors. 
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Figure  4.19:  Mach  0.35,  5000  ft.  Altitude,  Min  Def  CA  Mode 


Directional  command  variable  (deg/sec) 


Figure  4.20;  Mach  0.35,  5000  ft.  Altitude,  Min  Def  CA  Mode 
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Figure  4.21:  Mach  0.35,  5000  ft.  Altitude,  Min  Def  CA  Mode 


Left  AMT  (deg) 


-  CMD 

—  DEF 
••'PREF 

Right  AMT  (deg) 


Symmetric  AMT  (deg) 


Asymmetric  AMT  (deg) 


Figure  4.22:  Mach  0.35,  5000  ft.  Altitude,  Min  Def  CA  Mode 
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Figure  4.25:  Mach  0.35,  5000  ft.  Altitude,  Min  Def  CA  Mode 
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Roll  acceleration  —  aero  model  error  (deg/(sec)^2) 


Roll  acceleration  -  actuator  model  error  (deg/(sec)'^2) 


Figure  4.27;  Mach  0.35,  5000  ft.  Altitude,  Min  Def  CA  Mode 


Pitch  acceleration  -  aero  model  error  (deg/(sec)^2) 


Pitch  acceleration  —  actuator  model  error  (deg/(sec)'^2) 


Figure  4.28:  Mach  0.35,  5000  ft.  Altitude,  Min  Def  CA  Mode 
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Yaw  acceleration  —  aero  model  error  (deg/(sec)'^2) 


Yaw  acceleration  —  actuator  model  error  (deg/(sec)^2) 


Figure  4.29:  Mach  0.35,  5000  ft.  Altitude,  Min  Def  CA  Mode 
Minimum  drag  mode 

The  minimum  drag  control  allocation  mode  time  histories  are  contained  in  Figs.(4.30)-(4.40). 
Note  that  there  is  more  coupling  between  command  variable  axes,  and  angle  of  slip  regulation 
is  good.  Elevons  and  yaw  thrust  vectoring  are  the  control  effectors  that  deviate  from  their 
preferences  the  most  with  all  other  surfaces  only  deviating  from  their  preferences  at  the 
onset  of  stick  commands  when  insufficient  control  power  exists  to  achieve  the  control  law 
acceleration  commands.  The  leading  edge  flaps  do  not  deviate  from  their  preference,  but 
note  the  preference  is  non-zero.  More  directional  oscillations  are  prevalent  in  the  yaw  thrust 
vector,  directional  command  variable  and  side  acceleration  responses.  This  again  is  due 
to  the  lower  bandwidth  actuator  of  the  thrust  vector  nozzle  compared  to  the  aerodynamic 
effectors.  There  is  almost  no  on-board  actuator  model  errors  and  only  slight  aerodynamic 
model  errors  in  the  lateral  and  directional  axes.  There  is  significantly  more  longitudinal 
aerodynamic  errors  due  to  the  unmodeled  interactions  between  the  leading  edge  flaps  and 
elevons  which  leads  to  the  command  variable  coupling. 
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Roll  acceleration  —  actuator  model  error  (deg/(sec)^2) 
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Figure  4.38:  Mach  0.35,  5000  ft.  Altitude,  Min  Drag  CA  Mode 


Pitch  acceleration  -  aero  model  error  (deg/(sec)^2) 


Pitch  acceleration  —  actuator  model  error  (deg/(sec)^2) 


Figure  4.39:  Mach  0.35,  5000  ft.  Altitude,  Min  Drag  CA  Mode 
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Yaw  acceleration  —  aero  model  error  (deg/(sec)'^2) 


Yaw  acceleration  —  actuator  model  error  (deg/(sec)^2) 


Figure  4.40:  Mach  0.35,  5000  ft.  Altitude,  Min  Drag  CA  Mode 
Minimum  wing  load  mode 

The  minimum  wing  load  control  allocation  mode  time  histories  are  contained  in  Figs.  (4.41)- 
(4.51).  Note  that  there  is  much  more  coupling  between  command  variable  axes,  and  angle 
of  slip  regulation  is  good.  The  spoiler  slot  deflectors  and  pitch  flaps  deviate  from  their 
preferences  more  in  this  mode  to  make  up  for  stronger  preference  adherence  of  the  all  moving 
tips  and  leading  edge  flaps.  All  moving  tips  and  leading  edge  flaps  only  deviate  from  their 
preferences  at  the  onset  of  stick  commands  when  insufficient  control  power  exists  to  achieve 
the  control  law  acceleration  commands  since  these  surfaces  produce  large  wing  loads.  Even 
more  directional  oscillations  prevalent  in  the  yaw  thrust  vector,  directional  command  variable 
and  side  acceleration  responses  exist.  This  is  due  to  the  lower  bandwidth  actuator  of  the 
thrust  vector  nozzle  compared  to  the  aerodynamic  effectors.  There  is  almost  no  on-board 
actuator  model  errors  and  only  slight  aerodynamic  model  errors  in  the  longitiidinal  and 
directional  axes.  More  lateral  aerodynamic  errors  are  prevalent  due  to  unmodeled  spoiler 
slot  deflector  and  elevon  interactions. 
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Lateral  command  variable  (deg/sec) 


Longitudinal  command  variable  (deg/sec) 


Figure  4.41:  Mach  0.35,  5000  ft.  Altitude,  Min  Load  CA  Mode 


Directional  command  variable  (deg/sec) 


Figure  4.42:  Mach  0.35,  5000  ft.  Altitude,  Min  Load  CA  Mode 
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Figure  4.43:  Mach  0.35,  5000  ft.  Altitude,  Min  Load  CA  Mode 
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Figure  4.44:  Mach  0.35,  5000  ft.  Altitude,  Min  Load  CA  Mode 
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Left  SSD  (deg) 


Right  SSD  (deg) 
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Asymmetric  SSD  (deg) 


Figure  4.45:  Mach  0.35,  5000  ft.  Altittide,  Min  Load  CA  Mode 
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Figure  4.46:  Mach  0.35,  5000  ft.  Altitude,  Min  Load  CA  Mode 
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Figure  4.49:  Mach  0.35,  5000  ft.  Altitude,  Min  Load  CA  Mode 


Pitch  acceleration  -  aero  model  error  (deg/(sec)^2) 


Yaw  acceleration  -  aero  model  error  (deg/(sec)'^2) 


Yaw  acceleration  -  actuator  model  error  (deg/(sec)''2) 


Figure  4.51:  Mach  0.35,  5000  ft.  Altitude,  Min  Load  CA  Mode 

Minimum  radar  signature  mode 

The  minimum  radar  signature  control  allocation  mode  time  histories  are  contained  in  Figs.  (4.52)- 
(4.62).  Note  that  there  is  only  slight  coupling  between  command  variable  axes,  and  angle 
of  slip  regulation  is  very  good.  Elevons  and  thrust  vectoring  are  the  most  active  control 
effectors  with  all  other  surfaces  only  deviating  from  their  preferences  at  the  onset  of  stick 
commands  when  insufficient  control  power  exists  to  achieve  the  control  law  acceleration 
commands.  The  spoiler  slot  deflectors  and  all  moving  tips  follow  their  preferences  closely 
since  these  effectors  increase  radar  signature  more  than  other  effectors.  There  are  still  slight 
directional  oscillations  prevalent  in  the  yaw  thrust  vector,  directional  command  variable  and 
side  acceleration  responses.  This  is  due  to  the  lower  bandwidth  actuator  of  the  thrust  vector 
nozzle  compared  to  the  aerodynamic  effectors.  There  is  almost  no  on-board  actuator  model 
errors  and  only  slight  aerodynamic  model  errors. 
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Figure  4.52:  Mach  0.35,  5000  ft.  Altitude,  Min  Signature  CA  Mode 


Directional  command  variable  (deg/sec) 


Figure  4.53:  Mach  0.35,  5000  ft.  Altitude,  Min  Signature  CA  Mode 
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Figure  4.54:  Mach  0.35,  5000  ft.  Altitude,  Min  Signature  CA  Mode 


Left  AMT  (deg)  Right  AMT  (deg) 


Symmetric  AMT  (deg) 


Asymmetric  AMT  (deg) 


Figure  4.55;  Mach  0.35,  5000  ft.  Altitude,  Min  Signature  CA  Mode 
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Figure  4.58:  Mach  0,35,  5000  ft.  Altitude,  Min  Signature  CA  Mode 
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Figure  4.59:  Mach  0.35,  5000  ft.  Altitude,  Min  Signature  CA  Mode 
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Figure  4.60:  Mach  0.35,  5000  ft.  Altitude,  Min  Signature  CA  Mode 


Pitch  acceleration  —  aero  model  error  (deg/(sec)^2) 


Pitch  acceleration  —  actuator  model  error  (deg/(sec)^2) 


Figure  4.61:  Mach  0.35,  5000  ft.  Altitude,  Min  Signatiire  CA  Mode 


Air  Force  Research  Laboratory 


Yaw  acceleration  -  aero  model  error  (deg/(sec)'^2) 


Figure  4.62:  Mach  0.35,  5000  ft.  Altitxide,  Min  Signature  CA  Mode 

Injection  mode 

The  injection  control  allocation  mode  time  histories  are  contained  in  Figs.(4.63)-(4.73).  Note 
that  there  is  almost  no  coupling  between  command  variable  axes,  and  angle  of  slip  regulation 
is  very  good.  All  surfaces  are  active  and  follow  their  preferences  very  closely  since  the 
preferences  were  generated  with  knowledge  of  the  control  law  command.  There  are  no 
directional  oscillations.  There  is  almost  no  on-board  actuator  model  errors  and  only  slight 
aerodynamic  model  errors. 
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Lateral  command  variable  (deg/sec) 


Longitudinal  command  variable  (deg/sec) 


Figure  4.63:  Mach  0.35,  5000  ft.  Altitude,  Injection  CA  Mode 
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Figure  4.64:  Mach  0.35,  5000  ft.  Altitude,  Injection  CA  Mode 
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Figure  4.65:  Mach  0.35,  5000  ft.  Altitude,  Injection  CA  Mode 
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Figure  4.66:  Mach  0.35,  5000  ft.  Altitude,  Injection  CA  Mode 
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Figure  4.67:  Mach  0.35,  5000  ft.  Altitude,  Injection  CA  Mode 
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Figure  4.68:  Mach  0.35,  5000  ft.  Altitude,  Injection  CA  Mode 
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Yaw  acceleration  —  aero  model  error  (deg/(sec)'^2) 


Figure  4.73;  Mach  0.35,  5000  ft.  Altitude,  Injection  CA  Mode 


Minimum  thrust  vector  mode 

The  minimum  thrust  vector  control  allocation  mode  time  histories  are  contained  in  Figs. (4.74)- 
(4.84).  This  mode  was  a  late  addition  to  the  control  law  attempting  to  alleviate  the  direc¬ 
tional  oscillations  prevalent  in  the  other  modes.  Note  that  there  is  only  slight  coupling 
between  command  variable  axes,  and  angle  of  slip  regulation  is  very  good.  Thrust  vectoring 
only  deviates  from  the  preferences  at  the  onset  of  stick  commands  when  insufficient  control 
power  exists  to  achieve  the  control  law  acceleration  commands.  There  are  no  directional 
oscillations  which  led  to  the  nominal  choice  of  this  mode  at  the  low  dynamic  pressure  flight 
condition.  There  is  almost  no  on-board  actuator  model  errors  and  only  slight  aerodynamic 
model  errors. 
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Lateral  command  variable  (deg/sec) 


Longitudinal  command  variable  (deg/sec) 


Figure  4.74:  Mach  0.35,  5000  ft.  Altitude,  Min  Vec  CA  Mode 


Directional  command  variable  (deg/sec) 


Figure  4.75:  Mach  0.35,  5000  ft.  Altitude,  Min  Vec  CA  Mode 
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Figure  4.76:  Mach  0.35,  5000  ft.  Altitude,  Min  Vec  CA  Mode 
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Figure  4.77:  Mach  0.35,  5000  ft.  Altitude,  Min  Vec  CA  Mode 
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Roll  acceleration  -  actuator  model  error  (deg/(sec)^2) 


Figure  4.82:  Mach  0.35,  5000  ft.  Altitude,  Min  Vec  CA  Mode 


Pitch  acceleration  -  aero  model  error  (deg/(sec)'^2) 


Pitch  acceleration  -  actuator  model  error  (deg/(sec)''2) 


Figure  4.83:  Mach  0.35,  5000  ft.  Altitude,  Min  Vec  CA  Mode 
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Yaw  acceleration  —  aero  model  error  (deg/(sec)''2) 


Yaw  acceleration  -  actuator  model  error  (deg/(sec)''2) 


Figure  4.84;  Mach  0.35,  5000  ft.  Altitude,  Min  Vec  CA  Mode 


Summary 

The  various  control  allocation  modes  minimize  their  corresponding  objectives.  The  linear 
program  control  allocation  deviates  the  minimum  number  of  control  effectors  from  their 
preferences  to  achieve  the  control  law  command.  The  injection  mode  results  in  almost  no 
residual  aircraft  motion  and  decoupled  command  variable  response.  Directional  oscillations 
occur  for  those  modes  that  use  thrust  vectoring  as  the  primary  yaw  control  effector.  This  is 
attributed  to  the  fact  that  the  thrust  vectoring  nozzle  has  the  lowest  bandwidth  actuator  of 
all  the  actuators.  The  low  bandwidth  vector  actuator  reduces  the  stability  margins  at  this 
flight  condition.  Further  recall  that  the  zero  dynamics  were  least  stable  and  most  sensitive 
to  control  allocation  at  this  flight  condition.  The  nominal  mode  for  this  flight  condition  was 
chosen  to  be  minimum  vectoring  due  to  these  factors.  Command  variable  response  coupling 
was  prevalent,  in  those  modes  that  led  to  large  leading  edge  flap  and  spoiler  slot  deflector 
deflections,  due  to  unmodeled  interactions  of  these  surfaces  with  the  elevens. 

4.4.2  Final  control  law  analysis 

This  section  contains  batch  simiilation  responses  of  loaded  roll  and  full  stick  roll  maneiivers 
with  the  control  laws  in  the  final  configuration  at  each  of  the  four  analysis  flight  conditions. 

Low  dynamic  pressure  flight  condition,  Mach  0.35  5,000  ft.  altitude 

This  section  contains  simulation  responses  at  the  low  dynamic  pressure  flight  condition.  The 
control  allocation  mode  at  this  condition  minimizes  thrust  vectoring. 
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The  responses  in  Figs.(4.85)-(4.95)  are  for  the  loaded  roll  manenver.  These  responses  are 
identical  to  the  responses  in  Figs.(4.74)-(4.84).  There  is  only  slight  coupling  between  com¬ 
mand  variable  axes,  and  angle  of  slip  regulation  is  very  good.  Thrust  vectoring  only  deviates 
from  the  preferences  at  the  onset  of  stick  commands  when  insufficient  control  power  exists 
to  achieve  the  control  law  acceleration  commands.  There  is  almost  no  on-board  actuator 
model  errors  and  only  slight  aerodynamic  model  errors. 
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Longitudinal  command  variable  (deg/sec) 


Figure  4.85:  Mach  0.35,  5,000  ft.  Altitude,  A4in  Thrust  Vector  CA  Mode,  Loaded  Roll 


Directional  command  variable  (deg/sec) 


Figure  4.86:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Loaded  Roll 
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Figure  4.87:  Mach  0.35.  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Loaded  Roll 


Left  AMT  (deg)  Right  AMT  (deg) 


Figure  4.88:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Loaded  Roll 
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Left  SSD  (deg)  Right  SSD  (deg) 


Figure  4.89;  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Loaded  Roll 


Left  OBLF  (deg)  Right  OBLF  (deg) 


Figure  4.90:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Loaded  Roll 
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Figure  4.91:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Loaded  Roll 
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Figure  4.92:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Loaded  Roll 


Accelerations  (g) 
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Roll  acceleration  -  aero  model  error  (deg/(sec)''2) 


Roll  acceleration  —  actuator  model  error  (deg/(sec)'^2) 


Figure  4.93:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Loaded  Roll 


Figure  4.94:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Loaded  Roll 
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Figure  4.95:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Loaded  Roll 

The  responses  in  Figs. (4.96)-(4. 106)  are  for  the  full  stick  roll  maneuver.  The  lateral 
command  variable  tracks  the  command  well  and  slight  directional  coupling  exists.  Angle 
of  slip  regulation  is  very  good.  Longitudinal  coupling  caused  by  aerodynamic  model  errors 
degrades  longitudinal  command  variable  regulation.  The  aerodynamic  model  errors  are  due 
to  unmodeled  spoiler  slot  deflector  and  leading  edge  flap  interactions  with  the  elevons.  Thrust 
vector  deflections  only  deviate  from  their  preferences  at  the  onset  of  stick  commands  when 
insufficient  control  power  exists  to  achieve  the  control  law  acceleration  commands.  There 
are  roll  acceleration  errors  due  to  unmodeled  actuator  dynamics. 


Yaw  acceleration  —  aero  model  error  (deg/(sec)''2) 


Yaw  acceleration  -  actuator  model  error  (deg/(sec)''2) 
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Lateral  command  variable  (deg/sec) 


Figure  4.96:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Full  Stick  Roll 


Directional  command  variable  (deg/sec) 


Figure  4.97:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Full  Stick  Roll 
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Left  ELEV  (deg)  Right  ELEV  (deg) 


Figure  4.98:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  F\ill  Stick  Roll 


Left  AMT  (deg)  Right  AMT  (deg) 


Figure  4.99:  Mach  0.35,  5,000  ft.  Altitiide,  Min  Thrust  Vector  CA  Mode,  Full  Stick  Roll 
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Figure  4.100:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Full  Stick  Roll 
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Figure  4.101:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Full  Stick  Roll 
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Figure  4.102:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Full  Stick  Roll 
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Figure  4.103:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Full  Stick  Roll 
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Roll  acceleration  -  aero  model  error  (deg/<sec)'"2) 


Roll  acceleration  ~  actuator  model  error  (deg/(sec)^2) 


Figure  4.104:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Full  Stick  Roll 


Pitch  acceleration  -  aero  model  error  <deg/(sec)^2) 


Figure  4.105:  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Full  Stick  Roll 
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Yaw  acceleration  —  aero  moctel  error  (deg/(sec)'^2) 


Figure  4.106;  Mach  0.35,  5,000  ft.  Altitude,  Min  Thrust  Vector  CA  Mode,  Full  Stick  Roll 

Subsonic  ingress  flight  condition,  Mach  0.90  1,000  ft.  altitude 

This  section  contains  simulation  responses  at  the  subsonic  ingress  flight  condition.  The 
control  allocation  mode  at  this  condition  minimizes  radar  signature. 

The  responses  in  Figs.(4.107)-(4.117)  are  for  the  loaded  roll  maneuver.  There  is  only 
slight  coupling  between  command  variable  axes,  and  angle  of  slip  regulation  is  very  good. 
The  all  moving  tips,  spoiler  slot  deflectors,  leading  edge  flaps  and  pitch  flap  are  not  used 
for  this  maneuver  since  there  is  plenty  of  control  power,  and  use  of  these  effectors  increases 
radar  signature.  There  is  no  on-board  actuator  model  errors  and  only  slight  aerodynamic 
model  errors. 
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Figure  4.107:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Loaded  Roll 


Directional  command  variable  (deg/sec) 


-  CMD 
-  NCV 


Aero  angles  (deg) 


- AOA 

- AOS 


Figure  4.108:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Loaded  Roll 
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Figure  4.109;  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Loaded  Roll 


Figure  4.110:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Loaded  Roll 
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Left  SSD  (deg)  Right  SSD  (deg) 


Figure  4.111:  Mach  0.90,  1000  ft!  Altitude,  Min  Signature  CA  Mode,  Loaded  Roll 
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Figure  4.112:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Loaded  Roll 
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Figure  4.114:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Loaded  Roll 
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Roll  acceleration  -  aero  model  error  (deg/(sec)'^2) 


Figure  4.115:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Loaded  Roll 


Pitch  acceleration  -  aero  model  error  (deg/(sec)''2) 


Figure  4.116:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Loaded  Roll 
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Yaw  acceleration  -  aero  model  error  (deg/(sec)''2) 


Figure  4.117:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Loaded  Roll 

The  responses  in  Figs.(4.118)-(4.128)  are  for  the  full  stick  roll  maneuver.  The  lateral 
command  variable  tracks  the  command  well,  and  angle  of  slip  regulation  is  good.  Longi¬ 
tudinal  and  directional  coupling  caused  by  aerodynamic  model  errors  degrades  command 
variable  regulation.  The  aerodynamic  model  errors  are  due  to  low  model  fidelity  of  the  con¬ 
trol  derivatives  at  high  speeds.  Recall  that  the  on-board  high  speed  control  derivative  model 
is  constructed  from  linear  interpolated  linear  models.  The  large  elevon  deflections  violate 
the  small  angle  assumptions  of  the  on-board  linear  models.  A  higher  fidelity  high  speed 
control  derivative  model  would  alleviate  the  aerodynamic  model  errors.  Spoiler  slot  deflec¬ 
tors  and  all  moving  tips  only  deviate  from  their  preferences  at  the  onset  of  stick  commands 
when  insufficient  control  power  exists  to  achieve  the  control  law  acceleration  commands  since 
these  effectors  increase  radar  signature.  There  are  almost  no  roll  acceleration  errors  due  to 
unmodeled  actuator  dynamics. 


Air  Force  Research  Laboratory 


116 


Lateral  command  variable  (deg/sec) 


Figure  4.118:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Full  Stick  Roll 


Directional  command  variable  (deg/sec) 


Figure  4.119:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Full  Stick  Roll 
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Figure  4.120:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Full  Stick  Roll 
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Figure  4.121:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Full  Stick  Roll 
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Figure  4.122:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Full  Stick  Roil 
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Figure  4.123:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Full  Stick  Roll 


Air  Force  Research  Laboratory 


119 


5 


PTV  (deg) 


5 


YTV  (deg) 


PFLAP  (deg) 


1 

■ 

j _ 

0 

w  r 

-1 

■ 

I 

-2 

J 

1 

1 

1 

1 

-3 

-4 

1 

Figure  4.124:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Full  Stick  Roll 


Figure  4.125:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Full  Stick  Roll 
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Figure  4.126:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Full  Stick  Roll 
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Figure  4.127:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Full  Stick  Roll 
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Figure  4.128:  Mach  0.90,  1000  ft.  Altitude,  Min  Signature  CA  Mode,  Full  Stick  Roll 

Supersonic  cruise  flight  condition,  Mach  1.20  22,500  ft.  altitude 

This  section  contains  simulation  responses  at  the  supersonic  cruise  flight  condition.  The 
control  allocation  mode  at  this  condition  minimizes  drag  which  is  modeled  the  same  as 
minimum  deflection  at  this  flight  condition. 

The  responses  in  Figs.(4.129)-(4.139)  are  for  the  loaded  roll  maneuver.  There  is  only 
slight  coupling  between  command  variable  axes,  and  angle  of  slip  regulation  is  very  good.  The 
elevens  and  yaw  thrust  vectoring  are  the  primary  control  effectors  used  for  this  maneuver. 
The  all  moving  tips,  spoiler  slot  deflectors,  pitch  thrust  vectoring  and  pitch  flap  are  only 
used  at  the  onset  of  stick  commands  where  there  is  insufficient  control  power  to  achieve  the 
control  law  command.  The  leading  edge  flaps  are  not  used  for  this  maneuver.  There  is  no 
on-board  actuator  model  errors  and  only  slight  aerodynamic  model  errors  due  to  low  fidelity 
on-board  control  derivative  models  at  high  speeds. 
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Figure  4.129:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Loaded  Roll 
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Figure  4.130:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Loaded  Roll 
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Figure  4.131:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Loaded  Roll 


Left  AMT  (deg)  Right  AMT  (deg) 


Figure  4.132:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Loaded  Roll 
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Figure  4.133:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Loaded  Roll 
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Figure  4.134:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Loaded  Roll 
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Figure  4.137:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Loaded  Roll 
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Figure  4.138:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Loaded  Roll 
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Yaw  acceleration  -  aero  model  error  (deg/(sec)''2) 
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Figure  4.139:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Loaded  Roll 

The  responses  in  Figs.  (4. 140)- (4. 150)  are  for  the  full  stick  roll  maneuver.  The  lateral 
command  variable  tracks  the  command  well  ,and  angle  of  slip  regulation  is  degraded.  Lon¬ 
gitudinal  and  directional  coupling  caused  by  aerodynamic  model  errors  degrades  command 
variable  regulation.  The  aerodynamic  model  errors  are  due  to  low  model  fidelity  of  the  con¬ 
trol  derivatives  at  high  speeds.  There  are  almost  no  roll  acceleration  errors  due  to  unmodeled 
actuator  dynamics. 
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Figure  4.140:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Full  Stick  Roll 
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Figure  4.141:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Full  Stick  Roll 
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Figure  4.142;  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Full  Stick  Roll 
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Figure  4.143:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Full  Stick  Roll 
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Figure  4.145:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Full  Stick  Roll 
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Figure  4.147:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Full  Stick  Roll 
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Figure  4.148:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Full  Stick  Roll 
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Figure  4.149:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Full  Stick  Roll 
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Yaw  acceleration  -  actuator  model  error  (deg/(sec)'^2) 


Figure  4.150:  Mach  1.2,  22,500  ft.  Altitude,  Min  Drag  CA  Mode,  Full  Stick  Roll 

Envelope  center  flight  condition,  Mach  0.60  15,000  ft.  altitude 

This  section  contains  simulation  responses  at  the  envelope  center  flight  condition.  The 
control  allocation  mode  at  this  condition  minimizes  wing  loads. 

The  responses  in  Figs.(4.151)-(4.161)  are  for  the  loaded  roll  maneuver.  There  is  only 
slight  coupling  between  command  variable  axes,  and  angle  of  slip  regiilation  is  very  good.  The 
elevons  and  yaw  thrust  vectoring  are  the  primary  control  effectors  used  for  this  maneuver. 
The  all  moving  tips,  leading  edge  flaps  and  spoiler  slot  deflectors  are  only  used  at  the  onset  of 
stick  commands  where  there  is  insufficient  control  power  to  achieve  the  control  law  command 
since  these  effectors  caiise  high  wing  loads.  There  are  no  on-board  actuator  model  errors 
and  only  slight  aerodynamic  model  errors. 
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Figure  4.151:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Loaded  Roll 


Directional  command  variable  (deg/sec) 


Figure  4.152:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Loaded  Roll 
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Figure  4.153:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Loaded  Roll 
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Figure  4.154:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Loaded  Roll 
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Figure  4.155:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Loaded  Roll 
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Figure  4.159:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Loaded  Roll 
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Figure  4.160:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Loaded  Roll 
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Yaw  acceleration  —  aero  model  error  (deg/(sec)'^2) 


Yaw  acceleration  —  actuator  model  error  (deg/(sec)^2) 


Figure  4.161:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Loaded  Roll 

The  responses  in  Figs.(4.162)-(4.172)  are  for  the  full  stick  roll  maneuver.  The  lateral 
command  variable  tracks  the  command  well,  and  angle  of  slip  regulation  is  degraded.  Lon¬ 
gitudinal  and  directional  coupling  caused  by  aerodynamic  model  errors  degrades  command 
variable  regulation.  The  aerodynamic  model  errors  are  due  unmodeled  spoiler  slot  deflector 
and  leading  edge  flap  interactions  with  the  elevons.  There  are  roll  acceleration  errors  due  to 
unmodeled  acttiator  dynamics. 
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Figure  4.162:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Full  Stick  Roll 


Directional  command  variable  (deg/sec) 


Figure  4.163:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Full  Stick  Roll 
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Figure  4.164:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Full  Stick  Roll 
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Figure  4.165:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Full  Stick  Roll 
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Figure  4.166:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Full  Stick  Roll 
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Figure  4.167:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Full  Stick  Roll 
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Figure  4.168:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Full  Stick  Roll 
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Figure  4.169:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Fnll  Stick  Roll 


Air  Force  Research  Laboratory 


144 


Yaw  acceleration  —  aero  model  error  (deg/(sec)^2) 


Figure  4.172:  Mach  0.6,  15,000  ft.  Altitude,  Min  Load  CA  Mode,  Full  Stick  Roll 


4.5  Piloted  simulation 

Evaluation  of  the  control  laws  in  a  real-time  environment  have  begun  to  evaluate  pilot  in  the 
loop  properties  of  the  control  system.  The  results  from  the  piloted  simulations  will  appear 
in  a  future  document. 
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Chapter  5 
Conclusions 


A  modular  flight  control  architecttire  was  developed  that  is  highly  re-usable  and  allows  inser¬ 
tion  of  technology  advances  in  different  aspects  of  control  theory  such  as  control  allocation 
and  modeling.  Dynamic  inversion  proved  to  be  a  intuitive  flight  control  approach  that  re¬ 
sulted  in  excellent  nominal  stability  and  command  augmentation.  Application  of  dynamic 
inversion  to  a  tailless  fighter  required  only  minor  extensions  from  previous  efforts  in  command 
variable  definitions.  The  angle  of  slip  gain  in  the  yaw  command  variable  had  to  be  increased, 
due  to  the  tailless  configuration,  for  sufficient  angle  of  slip  attenuation.  Robust  stability  of 
dynamic  inversion  was  good,  but  robust  performance  was  slightly  deficient.  Nonlinear  on¬ 
board  models  were  required  to  achieve  acceptable  control  law  performance.  Refinements  to 
increase  on-board  airframe  model  fidelity  are  still  reqTiired  to  achieve  improved  control  law 
performance.  The  ICE  configuration  101-3  control  effector  suite  of  redundant,  multi-axis, 
nonlinear  control  effectors  was  challenging.  A  multi-branch  linear  program  based  control 
allocation  algorithm  was  developed  that  utilizes  all  available  (affine)  control  power  without 
violating  actuator  limits.  However,  the  control  affine  assumption  inherent  in  the  control 
allocation  problem  prohibits  directly  accounting  for  nonlinearities  such  as  control  effector 
interactions.  Spoiler  slot  deflector  and  elevon  interactions  were  problematic  and  solved  by 
limiting  spoiler  slot  deflector  deflections.  Actuator  bandwidth  limitations  are  not  accounted 
for  by  the  control  allocation  algorithm  due  to  the  static  nature  of  the  control  allocation 
problem  formulation.  Neglecting  actuator  bandwidths  became  problematic  due  to  inconsis¬ 
tent  actuator  bandwidths  between  control  effectors.  Use  of  thrust  vectoring  was  restricted 
at  low  dynamic  pressures  due  to  the  low  phase  crossover  frequency  of  its  actuator  compared 
to  the  aerodynamic  effector  actuators.  The  control  allocation  method  easily  integrated  mis¬ 
sion  segment  objectives,  such  as  minimum  drag,  and  optimized  these  objectives.  Eigenvalue 
and  simulation  analysis  showed  that  poorly  designed  control  allocation  methods  may  reduce 
stability  margins  or  even  destabilize  the  aircraft.  An  innovative  null  space  injection  method 
was  implemented  as  a  control  allocation  mode  to  assist  in  a  future  reconfiguration  capa¬ 
bility  through  adaptation.  The  injection  mode  was  very  successfvil  in  that  it  commanded 
decorrelated  effector  positions  without  degrading  aircraft  response. 
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